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ABSTRACT 


An experimental and theoretical study was made of a mono- 
pole antenna mounted on a finite ground plane located above an infinite 
ground. Circular flat discs and hemispheres were used for the finite 
ground planes. Experiments were conducted over a 6 to 1 frequency 
range, with the length of monopole antenna fixed at a quarter wave- 
length at the unper end of the frequency band. The radii of finite 
ground planes used for experiments were generally less than a quarter 
wavelength. 

Input impedances were measured as a function of frequency 
at the base of the antenna, using as variables the radius and the locations 
of the ground plane with respect to an infinite ground below. Scale 
models were used to obtain measurements of radiation patterns and 
antenna current distributions. Results of these measurements are 
presented graphically. 

A theoretical analysis was also made of a monopole with a 
hemispherical grovnd plane on an infinite ground. Far-zone electro- 
magnetic fields were calculated as a function of both the ground plane 
radius and the frequency. Radiation resistances were also calculated. 

The results of the study indicate that the radiation resistance 
of an electrically short antenna may be increased significantly by 


locating it on a small ground plane above the infinite ground, rather 


than directly on the infinite ground. This conclusion opens the way for 


a more efficient utilization of receiving and transmitting antennas on 
ground- bised vehicles. The impedance characteristics of the antenna 


system are: such as to facilitate its operation over a broad frequency 


range. 


FOREWORD 


This report was prepared by the Cooley Electronics Labora- 
tory of The University of Michigan under United States Army Electronics 
Command Contract No. DA 28-043-AMC-02246(E), Project No. 5A6-79191- 
1902-02-24, "Improved Antenna Techniques Study." 

The research under this contract consists in part of an 
investigation to develop highly efficient remotely tuned impedance match- 
ing coupling networks for electrically short monopoles. 

The material reported herein represents a summary of a 
theoretical and experimental study which was made to determine the input 
impedance and radiation characteristics of an electrically short monopole 
over a small ground plane located at various distances above natural 


ground. 
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CHAPTER I 
INTRODUCTION 


Monopole and dipole antennas over an infinitely large con- 
ducting ground plane are two of the oldest types of antennas being usea. 
Their performance characteristics have been thoroughly studied in the 
past. However, the exact behavior of the antenna over a conducting 
ground plane of a finite size, and both antenna and finite ground plane 
over infinite ground with finite -onductivity has not been studied here- 


tofore. 


1.1 Statement of the Problem 

The purpose of this study is to investigate the properties 
of a vertical monopole antenna and a finite ground plane located above 
an infinite ground plane. Particular emphasis, in this investigation, 
is placed upon the study of antenna systems having circular ground planes 
ni which the dimensions of both the ground-plane diameters and the ground- 
plane heights above infinite ground are smaller than the wavelengths 
of interest. 

The investigation of the prepertics of this antenna system 
provic«s a basis for a detailed comparison with the conventional antenna 
syste! of a monopole antenna over an infinite ground plane and for de- 
termi: ation of the relative advantages and limitations of this system. 


In ad tition, the results of this study will be applied to an investigation 


of possible performance improvement for a monopole antenna mounted 


on a ground-based vehicle, since there is reason to believe that the an- 
tenna system proposed for study approximates the vehicular-mounted 


antenna. 


1.2 Topics of Investigation 


Since the purpose of this paper is to report a study of the 
vertical monopole with a finite ground plane in the presence of an infi- 
nite ground plane, the following topics are explored: 

(1) The experimental measurement of the input impedance 
of the antenna was carried out over the frequency range of 5 MHz to 
30 MHz (a 6 to 1 band) and as a funct ion both of the ground-plane size 
and of the ground-plane location, relative to an infinite ground. 

(2) The experimental measurement of current distribution 
on the antenna over the frequency range described above was carried 
out with the ground plane size and its location as the variables. 

(3) The radiation pattern was measured and its dependence 
on the ground plane size and location was studied. 

(4) A number of theoretical models was examined to deter- 
mine a suitable approximation to the given system. Results o.tained from 
the appropriate theoretical model was then c»mpared with the experimen- 
tal measurements. 


One of the theoretical models which was studied is a mono- 


pole placed on a semi-spherical conductor, as shown in Fig. 1. 


antenna 


conducting 
semi-sphere 


-> — fe 


infinite ground 


image 


Fig. 1.1. Theoretical models 


1.3 Review of the Literature 

Many studies of the properties of vertical monopole antennas 
mounted both on infinite and on finite ground planes have previously been 
reported in the literature. 

Bardeen (Ref. 1), in 1930, studied the diffraction of a cir- 
cularly symmetrical electromagnetic wave by a circular disc of infinite 
conductivity. The result of his study has been used in det ermining the 
power flow into the earth below a vertical antenna that is grounded by 
a circular disc lying on the ies: In 1945, Brown and Woodward 
(Ref. 2), in an experimental study, measured the resistance and reac- 
tance of a cylindrical antenna operated against ground. The ground im- 
mediately below the antenna had a buried metallic screen of a finite 
size whose top surface was flush with the rest of the natural ground. 

A particular emphasis was given to a study of impedance behavior as a 
function of the antenna length to diameter ratio. Terminal conditions 
such as capacitance of the base of the antenna to ground were considered. 

Later, Meier and Summers (Ref. 2) (1949) performed an 
experimental study of the impedance characteristics of vertical antennas 
mounted on finite ground planes. Leitner and Spence (Ref. 3) (1951) 
confirmed some of these experimental results through the theoretical 
study of a quarter-wavelength monopole on a finite, circular-disc ground 
plane. They assumed a sinusoidal current distribution to exist on the 
antenna. The results of this study showed a marked dependence ot the 


antenna radiation resistance upon the diameter of the disc employed. 


This dependence was particularly pronounced for small ground planes. 

In 1951, Storer (Ref. 4) used a variational method to readily obtain the 
expression for the dependence of the antenna impedance on the ground 
plane diameter. Tang (Ref. 5), in 1962, further studied the radiation 
pattern of this system by using the Babinet principle and the Wiener-Hopf 
technique. 

In each of the above papers, the monopole was considered to be 
above a finite ground plane suspended in a free space with no direct effect 
from the infinite ground plane. Solutions by Leitner and Spence (Ref. 3) 
are mainly for the ground plane diameter smaller than the wavelength 
of interest. On the other hand, Tang and Storer's solutions can best 
be applied to the ground plane larger than ten wavelengths in diameter. 

Wait and Pope (Ref. 6) in 1955, studied the input resistance 
of low-frequency, monopole antennas with the radial-wire earth system. 
The conducting wires were placed radially over the infinite natural ground 
and the changes in the input resistance AR from the input resistance of 
a monopole with an infinite conducting ground plane was calculated. 

Brown (Ref. 7), in 1937, and Abbott (Ref. 8), in 1952, carried 


out earlier studies of the radial wire system. 


1.4 Report Organization 
(1) Chapter 2 contains the development of the experimental 


procedures for the input impedance measurement of the combined antenna 


and for the ground plane. The resulcs of both actual and scale model 
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impedance measurements are included. Some discussion of the losses 
involved in the real part of the impedance is included. 

(2) Chapter 3 contains a discussion of the technique for 
making current measurements on the antennas. Results of the measure- 
ments are shown in graphic form. Also, a comparison is made with 
the sinusoidal current distribution that is ordinarily assumed for a mono- 
pole antenna over an infinite conducting ground plane. 

(3) Chapter 4 contains a discussion of the technique of ra- 
diation-pattern measurements. The pertinent assumptions made in ob- 
taining these patterns are discussed. The results of the experiment are 
displayed graphically in polar plots that show the dependence of side-lobe 
levels on the ground plane size and its location. 

(4) Chapter 5 contains an analysis of the theoretical model 
outlined in the topics of investigation. The far-field electric and mag- 
netic fields are calculated for a given current distribution. The radiation 
resistances are also calculated for different antennas and the results are 
compared with the experimental ones. The raciation patterns obtained 
from this theory are also compared with the experimental results. 

(5) Chapter 6 contains a summary and conclusions. Sug~- 
gestions for further research are presented and practical applications 
usi ng the result of this study are discussed. 

The measurement procedure is described in Section 2. 2 


in which the steps taken to ensure accuracy in the measurement of small 


input resistances are given. The results determined using this procedure 
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are reported in some detail in Section 2. 3, and compared with the input 
impedance of a monopole on infinite ground plane. Some additional 
measurements, which are set out in Section 2.4, were carried out with 
scale models (1) of the monopole over a small ground plane and both 
above the infinite ground and (2) of the monopole on a hemispherical 
ground plane above an infinite ground. Experimental measurements on 
the latter model were carried out and were compared with the theoretical 


results in Chapter 5. 


CHAPTER II 


INPUT IMPEDANCE MEASUREMENT 


2.1 Introduction 

This chapter contains a discussion of the experimental 
work carried out to measure the impedance of a monopole antenna 
that was a quarter-wavelength long at 30 MHz and that was located on 
a finite ground that, in turn, was located above natural ground. The 
impedance was measured as a function both of ground plane size and 
of ground plane location. 

Theretofore, no extensive study of such a configuration has 
been reported. The input impedance characteristics of an antenna 
which is mounted on an infinite ground plane has been reported by 
Brown and Woodward (Ref. 2) for a wide variety of antenna lengths 
and diameters. Both experimental and theoretical studies have been 
conducted on the input impedance of a monopole antenna mounted on a 
finite ground plane, having diameters ranging down to wavelength. 
These studies did not, however, consider the effect of an infinite 


grourd below the finite ground plane. 


2,2 Experimental Measurement Procedure 
Measurements were conducted with an antenna made of 
copper pipe One~quarter inch in diameter, which is a quarter-wave- 


length long at 30 MHz. The antenna was mounted on a non-conducting 


variable height test stand, shown in Fig. 2.1, With this stand, one 


Fig. 2.1. Antenna on the variable 
height ground plane 


16 


can take measurements at heights of from zero to 5 meters above 
ground within continuous changes, and at the ground plane diameters 
of from zero to 5 meters. 

In order to observe the impedance variation with size and 
location of the ground plane as measured in terms of wavelength, the 
size and the location of the finite ground plane was varied in four 
steps ranging from one-sixteenth of a wavelength to one-half wavelength 
at 30 MHz 

The block diagram in Fig. 2. 2 shows the test set-up. The 
test equipment was placed just below the ground plane, to minimize 
the reflection of electromagnetic waves from the equipment. The 
entire experiment was conducted in a flat field where the nearest 
building was at least 10 wavelengths away from the test set-up. (See 
Fig. 2. 3). 

Four sizes of ground plane were used in the experiments. 
A highly conducting 5 mil thick aluminum sheet was chosen for the 
ground plane material, after considering the skin depth at the frequen- 
cies where the measurements were taken. 

Because of the provision for variable height and diameter, 
the total weight of the ground plane is an important factor in designing 
the experimental set-up. The depth of penetration, 6 , for aluminum 
at frequencies of 5 MHz and 30 MHz is 1.45 mils and 0. 594 mil, re- 
spectively. The thickness chosen was the smallest,but was much 


thicker than the skin depth in aluminum sheet ground plane. A 5 mil 
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Antenna 


_ Finite Ground Plane 


External Capacitor 


General 


Collins 
HP 606A Radio 916A R-390/URR 
Signal RF Bridge Receiver 
Generator 


Fig. 2.2. Block diagram showing the impedance 
measurement set-up 
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Fig. 2.3. Antenna test site 
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Fig, 2.4. A 2.5 meter monopole antenna ona 2.5 
meter diameter ground plane supported 
by a styrofoam sheet 
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thick roll of aluminum sheeting is chosen for this purpose. 

To support the flexible ground plane, a hexagonal structure 
was built and it was covered with 4cm_ thick styrofoam sheet, which 
has a dielectric constant nearly equal to air.(Fig. 2.4). The ground 
plane was hoisted up to the proper height with ropes strung over rollers 
located at the top of the six supporting posts. 

Initial measurements were carried out using a length of 
50-ohm coaxial line to connect the base of the antenna to the :mpedance 
bridge. Readings taken with the bridge were referred to the antenna 
base by means of Smith charts. However, analysis shows appreciable 
inaccuracy to exist at low frequencies, where the antenna impedance Zr, is high. 
This inaccuracy was due to inaccuracies in the impedance transform- 
ation from the bridge terminals to the antenna base. In transmission 


line theory, an impedance Z d ata point d meters away from the 


load impedance a is 
Z. +jZ, tanwed 
L 0 r 
) ne (2. 1) 
. : Z,+jZ tan 22 d 
a ¢ x 


where Zy is the characteristic impedance of the line and A is the 
wavelength. If the load impedance is much greater thar te character- 
istic impedance a ae 25 at any frequency, then, this expression 


reduces to 
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= -jZ,cot-2= d (2.2) 


Z oot 


d 


which is independent of the load impedance Zr: Therefore, unless the 


characteristic impedance Zo is comparable in magnitude to Z a 


L’ 


small change in Z, is difficult to detect for large load impedances. 


L 
However, if d is chosen to be near zero then the changes in Zr, can 


be easily detected even if Zr, >> Zo: 


From equation 2.1 


Z + Zo tan” a d 
Bik gS Mesa oe Ae 
(2 + j Zr tan +— d) (2. 5) 
and therefore 
AZ A 
eeaUS, we pans 1 
24 Zs Z Z 


where 


In general, 
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Fig. 2.5. Test equipment arrangement for 
the impedance measurement 


: 
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Fig. 2.6. Location of the impedance measurement 
bridge relative to the antenna base 
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whenever Zo = Zr: 


However, from the Eq. 2.3 when 9 = O orwhen d— 0, 
d . L (2. 6) 


Equation 2.6 indicates that whenever the load impedance 


Zr, is much greater than the characteristic impedance of the line 


Zo the changes in Zr 


line length d = a where n_ is an integer. 


can be detected much easier by choosing the 


The input impedance of a monopole above an infinite ground 
plane and that resonates at a certain frequency fo is capacitive at 
frequencies smaller than fo: The magnitude of the capacitive react- 
ance is a rapidly varying function of frequency. 

Therefore, at the lower end of the frequency band, the 
technique of using a line extension to translate the impedance has a 
large error; and, for this reason, it was decided to modify the ground- 
plane support structure to accommodate the impedance bridge directly 
beneath the antenna base. This modification reduced the li length 
from the antenna to bridge to about 10 cm and thereby required a 
negligible impedance transformation. The set-up is shown in Figs. 
2.5 and 2.6. 

Measurements made using this set-up are hampered, how- 


ever, by the limited range of the bridge. It is capable of recording 


only reactances of magnitude less than 5,00C onms divided by the 
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Fig. 2.7. Antenna located at 9 meters above the 
natural ground with 2.5 meter diameter 
ground plane 
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measurement frequency in megahertz. The antenna capacitive react- 
ance at 5 MHz is greater than 1000 ohms. However, a high-Q capac- 
itor of approximately 36 u uf placed in shunt with the antenna across 
its base terminals reduced the overall reactance to the range covered 
by the bridge, and was, therefore, used. 

To obtain the input impedance of the antenna from the 
data obtained with the external shunt capacitor, the capacitor was itself 
carefully calibrated over the frequency range of 5 MHz to 30 MHz. 
Using this equipment, the capacitor impedance, with a nominal value 
of 36 1 uf, was measured at every 2.5 Miiz. The resistive component 
and the reactive component were carefully recorded at each frequency. 
These data were later used to obtain the actual input impedance of 
the antenna. The resistive component of the capacitor, over the 
frequency range of measurement, was always less than one-half of an 
ohm. In order to remove the errors of manual calculation, a computer 
vrogram was written in which the calibration data for the external 
shunt cavac‘i.or and the measurement data of the experiment were 
fed in as raw data and the series resistive and reactive components 
of the input impedance were calculated on the BM 7090 Computer. 

The calibration of the capacitor was done periodically 


during the period of this experiment. 


2.3 Scale Model Impedance Measurements 


Since Maxwell's equations are linear ones, an electro- 


magnetic structure that has certain properties at a given frequency 
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will have identical properties at another frequency nf, provided that 
all linear dimensions of the structure are scaled by the ratio 1/n. 
Therefore, an antenna design which operates over a certain range of 
frequencies can be made to operate over any other range of frequencies 
without additional re-design, if an exact scaling of dimensions is 
accomplished. 

Aside from permitting one to transfer design relationships, 
it is convenient to use a practical size scale model for radiation pat- 
tern studies. Scaling is a relatively simple matter for most types 
of antennas. However, electromagnetic properties must also ke scaled, 


as well as linear dimensions. The scaling factors are shown below. 


Table 2.1. Scale factrrs 


Both the length and frequency can be scaled easily, the 
conductivity cannot. However, the conductivity scaling is important 
only through losses, and, since these are small for most antennas, 
the inability to scale the conductivity exactly is not a serious problem. 

An exact scale model would retain the exact radiation 


patte.n and input impedance of the full scale antenna. However, it is 
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not always possible to scale the model exactly, particularly since 
the transmission lines and screw fastenings, etc. , cannot be scaled. 
Slight discrev..ncies in scaling will usually affect the impedance 
properties much more than they affect the radiation properties. The 
reason is that any discrepancies in scaling usually gives a rise in 
near-field effect which are attenuated out at far-field. 

The model used for the antenna and ground plane studies 
here is intended, therefore, mainly for radiation pattern studies. The 
scale mode’ .s designed to be 1/40 of the original antenna system in 
linear dimensions. This scaling factor allows the new antenna and 
ground plane system to operate between 200 MHz and 1200 MHz, 
replacing the original 5 MHz to 30 MHz range. 

As was pointed out previously, the principal difficulty in 
scaling lies with the conductivity. It is assumed here that the copper 
used for the scale model ground plane and the aluminum ground plane 
used for the full size antenna system both have infinitely large conduct- 
ivity, so that the no and o are both infinitely large and therefore, 
that losses are negligible. Also, the large ground plane covered with 
aluminum foil to simulate the natural ground at the higher frequencies 
is assumed to influence the radiation pattern only to a minor extent. 

Although the scale models are built mainly for radiation 
pattern studies, they will give qualitative information on the impedance 
as a function of frequency. The impedance measurements are, there- 


fore, carried out in the laboratory on the scale model between 200 
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MHz and 1200 MHz using a HP 805A slotted line and a HP 415D 
VSWR meter for frequencies above 500 MHz, anda HP 803A VHF 
impedance bridge with a HP 417A VHF detector for the frequencies 
below 500 MHz. A difficulty of measuring an impedance Zr, of an 
unknown ioad when Zr is much greater than the characteristic 
impedance of a transmission line Z» was pointed out in Section 2. 2. 
As one way to eliminate this difficulty, the line length d from the 
position of the load Zr to the measuring point Z q Was shortened. 
However, on the scale model simulating an infinitely large natural 
ground, the measurement equipment could not be mounted right at 
the input terminal of the antenna, thereby vs."ucing d to zero. 

A 50-ohm solid copper outer conductor coaxial line of 
approximately 20 centimeters was Sed to connect the antenna to the 
measuring equipment. Figure 2.8 shows a schematic diagram of the 
experimental arrangement. Each measurement was performed. 
first with a short at the position of the load and then with the load 
connected. This technique accounts for the length of line used between 
the measurement equipment and the load, when a Smith chart is used 
to obtain the actual input impedance of an antenna. The measurenient 
performed in this manner is reliable at higher frequencies where 
2 = Zo than it is at the lower frequencies where Zr PP Zo: 


The results are shown in Figs. 2.15 and 2. 16. 


In a further study, another set of scale models were 


built using hemispherical ground planes, instead of flat circular 
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Fig. 2.8. Block diagram showing an impedance measure- 
ment set-up for a scale model antenna system 
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ground planes, to examine any similarities in the radiation patterns 

of the two systems. The frequencies selected were such that the 
diameter of the hemispnerical ground planes was less than a wave- 
length at 1476 MHz. This is, in effect, a 49.2:1 scale model for 
Simulating operation between 5 MHz and 30 MHz. Impedance measure- 
ments were also made on this model, as well as the radiation pattern 
studies reported in Section 4. 4. 

In studying the impedance of the scale models, the flat 
circular disc ground planes were placed at heights ay and ao above 
the simulated natural ground. The displacement ay is the radius 
of the semi-spherical ground plane. The displacement ay is the 
radius of the flat disc equal in surface area to that of the spherical 
ground plane. 

Figure 2.9 shows, in detail, the experimental arrange- 
ments of scale models. In the first case, Fig. 2.9(a), the displace- 
ment oi a finite ground plane above an infinite ground is the same as 
the radius of a hemispherical ground plane a. In the second case, 
Fig. 2. 9(b),the displacement, a , of a finite ground plane above an 
infinite ground is chosen to be equal to 0. 707 times the radius of a 
circular disc. With this value, the surface area of the hemisphere, 
27a” , is equal to the surface area of the circular disc. That is, 
in the first case, displacements above an infinite ground plane are 


equal, and in the second case the surface areas of the small ground 


equal. The object of this study was to determine the relative 
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impedance characteristics of the two antenna systems and to determine 
the dependence of the system upon ground plane surface area and upon 
its location with respect to the larger ground plane. The results of 

the laboratory impedance measurement on both scale models are shown 


in Figs. 2.15 and 2. 16. 


2.4 Results of Measurement 

Antenna input impedance measurements taken with 5, 2.5, 
1.25, and 0.625 meter diameter ground plane at heights.5, 2.5, 1.25, 
and 0.625 meter above natural ground are shown as a function of 
frequency in Figs. 2.11 to 2.14. The data are arranged to show the 
effect of ground plane sizes with their location fixed at constant level 
and the effect of ground plane locations with the size fixed. 

The capacitive reactance components are, in general, 
less dependent upon the ground plane height than upon its diameter. 
However, the resistive component varies with the height markedly 
and shows a peak at particular frequencies. The input impedance of 
a monopole antenna that is located over an infinite ground plane is 
shown in Figs. 2. 10(a) and 2. 10(b). These curves are derived from 
Brown and Woodward's (Ref. 2) experimental data where the length 
to diameter ratio of approximately 400 is used. Comparison between 
results of infinite ground planes and finite ground planes above natural 
ground shows the two are similiar in reactance but not in resistance. 

The results from measurements using actual size of 


antenna show that, for a fixed diameter disc ground plane, a resistive 
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Fig.2,9. Scale models used for impedance measurements 


Resistive Component (Chms) 


28 


60 ->—- 


5 10 15 20 25 30 
Frequency (MHz) 


Fig. 2. 10(a).Input resistence as a function of frequency 
for a quarter wavelength monopole at 30 MHz 
on an infinite ground plane 
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Fig. 2.10(b). Input reactance as a function of frequency 
for a quarter wavelength monopole at 30 MHz 
On an infinite ground plane. 
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peak occurs at higher frequencies as the location of this ground plane 
gets higher with respect to an infinite ground. Also, the peak values 
of the resistive components are decreasing as the fixed diameter size 
is increased. A general conclusion which can be drawn from these 
results is that the resonant peak occurs at a frequency where 
K(- +a) is constant. Therefore, as a diameter of the ground plane 
and a displacement above infinite ground plane gets larger, a peaking 
occurs at lower frequencies. As for the reactive components of the 
impedance, the reactance versus frequency curve is similar to that 
of a monopole above an infinite ground plane when the diameter of a 
disc is 5 meters. The displacement above an infinite ground does not 
seem to affect the reactive component as much as the diameter of a 
finite ground plane. However, evenfor D = 5 meters, reactance 
values are lower than that of a monopole antenna with the same length 
to diameter ratio on an infinite ground. 

A limited number of impedance measurements has been 
conducted using a scale model which is mainly constructed to take a 
far-zone radiation pattern of the antenna and ground plane system. 
The input terminals as well as infinite natural ground simulations 
have not been able to scale ideally. These discrepancies usually 
affect impedance measurements more than radiation patterns cbtained 
at far-zone area. Under these circumstances, only a qualitative 
comparison of the impedance has been possible. 


First, the actual model and the scale model antenna and 


ground plane system at a ground plane diameter equal to a half wave- 
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length at 30 MHz compared favorably in both real part and imaginary 
part. 

Secondly, the scale model antenna with a disc ground plane 
and a hemispherical ground plane comparison shows that neither the 
diameter of the two t,'pes of ground planes or the area of the ground 
planes has a basis of similarity in impedances. However, it is 
deduced qualitatively from the limited number of experimental data 
that the distance between the base of the antenna to the natural ground 
and the conductivity along this path has a greater bearing upon the 
resonance phenomena observed in the impedance measurements. In 
other words, k(a + 2) = const for a disc ground plane and ka = 
const for a semi~spherical ground plane will determine the resonance 
conditions. These constants seem duferent in general for different 
ratio of a to D. 

In addition to these observations, the determination of 
whether resonance peak observed in input resistance measurement is 
largely due to increase in radiation resistance or not will be shown 
in the rest of this chapter and in Chapter 5. 

In studying impedance measurement data, it is also 
observed that both real and imaginary part of the impedance behave, 
in most cases, such as in Fig. 2.17(c). Frequencies where peaking 
effect occurs are, of course, a function of a and D. Figure 2. 17(a) 


shows an input impedance as a function of frequency for a monopole 


on an infinite conducting ground plane and its equivalent circuit. 


Pe TT eee 


Fig. 2.17. 
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(a) 


An equivalent circuit of a monopole on a finite 
disc ground plane above an infinite ground 


Freq. 


Freq. 


Freq. 
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Figure 2. 17(b) shows a parallel resonance circuit whose impedance 
characteristics add up with Fig. 217(a) to give an impedance char- 
acteristic which is observed in the measurements. 

Therefore, it is possible to synthesize an equivalent 
circuit of a type shown in Fig. 2.17(c) to further study the effect of 
a finite ground plane upon an equivalent circuit. In this way, a cor- 


relation may be obtained between a, D and the circuit parameters. 


2.5 Copper Losses Due to the Antenna and the Ground Plane 


The real part of the antenna input impedance contains a 


part 


1 


R =—> 
re ye 


§ 5 Re (E x H*)- dS (2. 9) 
that is directly proportional to the radiated power of an antenna, fora 
constant input current. Consequently, it is important to separate 

the resistive component of the input impedance measurement data 

into the resistive loss, which is due to several causes, as explained 

in the following section, and the radiation resistance. This permits 
determination of whether the unusual variation of the input resistance 
as a function of frequency, which has been found in the input impedance 
measurements, is due to an increase in loss or to an increase in the 
radiation resistances at particular frequencies. 


The total antenna resistance is the sum of the several 


separate components 
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(1) Radiation resistance R. 

(2) Ground terminal resistance R 
(3) Resistance of tuning units R, 
(4) Resistance of equivalent insulation loss R. 

(5) Resistance of equivalent conductor losses R., 

(6) Transmission line loss Rn ; 

Among the six separate contributions, R, and Rk. are 
not present in this case because there are no tuning units attached when 
the measurements are taken and, as a low potential receiving antenna, 
no appreciable insulation losses are involved. 

The loss due to the transmission line (Ry) was evaded 
by ineasuring the impedance at the base directly below the ground 
plane, for the actual size and by taking a short circuit measurement 
with the shcrt placed at the input terminal of the scale model. The 
loss in the line was subtracted from the measured data. 

In the out-door measurement with the natural ground 
below the finite aluminum ground plane R, cannot be exactly calcu- 
lated without precise information of the conductivity and other para- 
meters of the dirt ground. The impedance measurement on the scale 
model when the natural ground is simulated by the aluminum foil 
enables both R, and R to be computed. 

The following analysis is to permit the calculation of R 
for both the actual model and the scale model and Ry of ihe scale 


model. 
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2.9.1 The Internal Impedance of the Plane Conductor. 


The current density resulting from the movement of charges in a 


conductor is given by Ohm's law: 
J = oF (2. 8) 


The constant o is the conductivity of the conductor, and the Maxwell's 


equation is 
Vxii = ob + 2 (2. 9) 


For a harmonically oscillating field with e! sa time dependence 
VxH = (o+ jue)E (2. 10) 


In the absence of free charges p, V° D = 0. Also, for most con- 
ductors, the displacement current dD /dt is negligibly small com- 
pared to the conduction current. 


Then, 


VxVxE = V(V: E)-V°E = vx (=28) - < uvxi 


(2, 11) 


vV*E = po 2 = jw pok 


One 
u 
Q 
rl 


Similarly, using J 


seal | 
iu 


jwuod (2. 12) 
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Conducting Surface 


Fig. 2.18. Plane solid conductor 


ee 
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For a plane conductor of infinite depth, with no field vari- 


ations along the length or width as shown iu Fig. 2. 18. 


ad i, ack 
= = jwyuoi, = T i 2. 13) 
where 
T = jwuo 
Since 
Vi = {i+ j) 
NE} 
and 


T = (1+ j) Vito (2. 14) 


where by definition 5 = set , is called depth of penetration 


Viipo 


of the field, or the skin depth. 


The solution of differential equation 2.13 is then given as 


aS, 
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-(1 +i), 
i,=zipe © (2. 15) 


The total current flowing in the plane conductor is found by integrating 


the current density i, from the surface to the infinite depth. Fora 


unit width, 
7) we) = Cr) ig 6 
be = ne ack = 75 (2. 16) 
‘0 
The electric field at the surface is Evo oo Therefore, the 


internal impedance per unit length and unit width is 


E 14j 
Be ee AO a (2. 17) 


If Z_ isdefinedas Z_ = R_+jwL., then 
s Ss Ss i 


Rg -+.f me (2. 18) 
s a z0) o 


R, is the resistance of the plane conductor for a unit 
length and unit width. For a finite area of conductor, the resistance is 
obtained by multiplying R. by length, and dividing by the width. 

For a circular ground piane with radial current distribution, 
the total surface resistance is obtained by muitiplying R, by the 


radius and dividing by the mean circumference of the plate. 


For an aluminum ground plane which has material constants 
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3.72 x 10° mhos/ meter 


CG = 
-7 P 
pe = 4nx10  henries/meter 
6 = 0.0826 meters 
vf 
{ = frequency in Herts 


the surface resistivity is computed to be 


R, = 3.26x 10°! VF (2. 19) 
which becomes 
at { = 30 MHz R, = 3.26x10°%3x10' = 1.79% 10°° ohms/ 
square 
= 5MUz = 3.26x107%5x 10° = 7.3x10% " 


When the aluminum ground plane has a diameter equal toa 
quarter-wavelength at 30 MHz, the radius r is 2.5 meters and the 


mean circumference is mr. The total surface resistance is, therefore, 


R. 
R = R eke Gee (2. 20) 


The numerical values computed at 5 MHz and 30 MHz 


become 
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-3 
. 79 - 
R = Leen ohms = 5.7.x 10 4 ohms at 30 MHz 
7 7.3x 1074 4 


=— S477 ohms = 2.32x 10 “ohms at 5MHz. 


The contribution of the ground plane surface resistance toward the 


input resistance is, therefore, negligible. 


2.5.2 Internal Impedance of a Conductor with a Circular 
Cross Section 
2.5.2.1 Current in a Wire of a Circular Cross Section. 
Let the current flowing on the antenna of a circular cross section, 
1/4 inch in diameter be assumed to flow mainly on the axial direction; 


rd AN : : : Be Bact 2 
1 = i 2. Also, no axial or circumferential variation is assumed. 


It 


Antenna 


Fig. 2.19. Current in a cylindrical wire 
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Then V’ i = jwuoi becomes in the cylindrical coordinate system, 


i di 
gk +p Gh = iousi, en 
If welet T* = -jwyuo 
d?i di 
a sis Ti, = 0 (2, 22) 
Z 


For a solid wire, the solution must be finite at r = 0. 


Therefore, it takes the form of 


ae AJ (Tr) (2, 23) 


Let 
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then 
ty (2. 25) 
— Jo(Tr9) ) 
Jol) 2. 26 
z 0 Jo(Try) et) 
Since 
T= -jwpo 
(2, 27) 
= Fone 22¢ MnO [wuo (1 - j) 
T jwuc = F = 5 
Since 
ee ae t 4 Raj 2.28 
“oi) Ber (v) + j Bei (v) (2. 28) 
J 
Where 


Vi 


Ber (v) = Real part of Jo ( . ) 
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Bei (v) = Imaginary part of Jof v 


Wr 


The current density in the axial direction can be written 


as now 


i) 


Wor \ (lar 
Ber( ) + j Bei 5 (2. 29) 


If the ratio of r/6 is large, the li,/ ig| plot will agree 


Closely with the plane conductor derivation of 


e Tat (2. 30) 


where Ty -r replaced x for the case of a plane conductor. 


Also, 


@H-dl =I and 2QerH = 1 (2. 31) 
0 ¢ r=rg 


From the Maxwell's equations, 


VxE = -jwuH (2. 32) 


and for the round wire 
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Since 


1 
goes oe Ot. = 2.34 
E, “ ¢  ¢ Ip(Tr) ( ) 


H¢ = - TT Tir) (2. 35) 


Consequently, 


Pack amr Jo (Tr) 


F iy TT TF) (2. 36) 
2.5.2.2. The Internal Impedance of a Round Wire. 
The internal impedance Z. is defined as 
T J,(Tr,) 
alien aie os ee) 
i mr god 9 (Tr, 


Using the formula that 


Ber (v) + j Bei(v) = (=) 


vi 
and 


(Fr) 
Ber'(v) + j Bei'(v) = aa [Ber (v) + j Bei wv) = —\4) 


dv 


tts eas bahec cashes! “ech. ena 
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The internal impedance becomes 


wid noah eat 


Z. = R+ jwL epee See een et (2. 39) 
: : ; ——— : 
i i var, Ber' q+ j Bei'q 


where 
swede. te 
R, aaa: = rf uo (2. 40) 
and 
Vir 
_ 0 
q= = (2. 41) 
or 
eae R, | Ber q Bei’ q - Beiq Ber'q | ohms/meter 
m9 (Ber' q)* + (Bei q)* 
(2, 42) 
pa a Ber q Ber' q + Bei q Bei' q ohms/ meter 
1 Pm | Ber’ g)* + (Bei a)” 


Using the same analysis except that the wall thickness is 
small enough with respect to the radius of the tube to be able to consider 
the tublar conductor as a flat conductor of finite thickness, the internal 


impedance can be found to be 


2 = sar _Cosh Td_ 


(2. 43) 


Ss Sinh td 
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where 


(1 + j) 
0 


d = thickness of the tubular wall 


Ther 


ohm/unit length 
(2. 44) 


ohm/ unit length 


(2, 45) 


R (2a/ ) (2a/ ,) 
__S_ | Sinh o/ + Sin 7 / | ohm/ unit length 


dar a: 6) - Cos| } 6) 


0 Cosh 


The tubular antenna used for impedance measurement has 


the following dimensions: 


radius = 3.17 x 10° meters 


wall thickness = 6.3 x 1074 meters 
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length of antenna = 2.5 meters 
The skin depth 5 computed for the copper at both ends of 


the frequency band are 


ll x 10° meters at 30 MHz 


Lo) 
u 


3x 10° meters at 5 MHz 


Therefore, the ratio 2d/65 for equation 2.45 becomes 


2d ‘6 11.4 at 30 MHz 


42.4 at 5MHz 


For large values of x,both sinh x andcosh x approach 
(1/ 2)e" and sinh x >. sin x andcosh x |. cos x. Equation 


2.45 then simplifies to 


/ 
Ros -Ss2e° wh. when aq) 5 > 1 (2. 46) 


For the copper tubular antenna 2.5 meters long the real 


part of the internal impedance becomes 


Ro 2 0.189 ohms at 30 MHz 


0.063 ohms at 5 MHz 


The contributton of the copper loss, as shown above, toward 


the antenna input impedance is also nep*vible. 
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{n this section, a numerical calculation of ground terminal 
resistance R, and conductor loss resistance Ro was made. In the 
actal antenna and ground plane system, the conductivity of the infinite 
ground bziow the finite ground must be known. However, calculations 
pe: formed here did ct wake the natural ground into account. The only 
ground considered is an aluminum ground plane used as a finite ground 
below the antenna. Contributions from the tubular antenna of 2.5 
meters in length, used for a monopole, to R, were also considered. 

In the frequency range of 5 to 30 MHz, cee the experiment was 


carried out, the ground loss from the largest ground plane (D = 5 


meters) used was computed to be less than 5.7 x 10° * ‘chines: “Phe tose 
from the tubular antenna conductor itself was less than 0. 2 ohms. | 
In the scale model, the ground plane used to simulate an 
infinite natural ground was a 4-foot diameter aluminum circular 
ground plane. The conductivity data and skin-depth data available for 
aluminum, both large and small ground plane loss were calculated. 
Also, a Conduction loss from an AWG No. 28 wire of 
6.25 cm in length used to simulate a monopole antenna was calculated. 
The overall loss due to R, and R for the scale model was well 
below 1 ohm and there was no obvious peaking effect due to R. and 
R as frequency was varied from 200 MHz to 1200 MHz. 
Therefore, it seems that the resonant peak observed in the 
real part of the antenna input impedance in the experiment originate 


from sources other than R, and Ro 
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CHAPTER II 


CURRENT MEASUREMENTS 


3.1 Introduction 


In this chapter, an extensive study of the current distribu- 
tion on a monopole antenna one quarter wavelength long at 30 MHz is 
reported as functions of a ground plane size and its lacation. King and 
Harrison (Ref. 9) computed analytically the current distribution of 
a symmetrically center-driven antenna. For this dipole, the current 
distribution was found to be extremely close to the sinusoidal distri- 
bution that is normally assumed, so long as the antenna length to di- 
ameter ratio is very large. However, the eifect of a small ground 
plane whose diameter is less than a half wavelength upon the antenna 
current when the system is placed above a” infinite ground has not 
been reported. 

The particular emphasis here is to show the differences 
between the current distribution on an antenna over an infinite ground 
plane and the same size of antenna over a finite ground displaced less 


than a wavelength above an infinite ground plane. 


3.2 Theory of Current Probe 


In order to detect the current densities on the antenna, a 
rectangular loop antenna whose dimensions are much smaller than a 
wavelength was used as a probe. The theory of operaticn of this loop 


as a probe has been studred by King and Whiteside (Ref.10 ). The 
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CHAPTER III 


CURRENT MEASUREMENTS 


3.1 Introduction 

In this chapter, an extensive study of the current distribu- 
tion on a monopole antenna one quarter wavelength long at 30 MHz is 
reported as functions of a ground plane size and its location. King and 
Harrison (Ref. 9) computed analytically the current distribution of 
a symmetrically center-driven antenna. For this dipole, the current 
distribution was found to be extremely close to the sinusoidal distri- 
bution that is normally assumed, so long as the antenna length to di- 
ameter ratio is very large. However, the eifect of a small ground 
plane whose diameter is less than a half wavelength upon the antenna 
current when the system is placed above a” infinite ground has not 
been reported. 

The particular emphasis here is to show the differences 
between the current distribution on an antenna over an infinite ground 
plane and the same size of antenna over a finite ground displaced less 


than a wavelength above an infinite ground plane. 


3.2 Theory of Current Probe 

In order to detect the current densities on the antenna, a 
rectangular loop antenna whose dimensions are much smaller than a 
wavelength was used as a probe. The theory of operaticn of this loop 


as a probe has been studied by King and Whiteside (Ref.10 ). The 
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components are also pres.vnt, even though small, which are proportional 


to the average electric field in the plane of the loop. 


When electro- magnetic fields are incident upon a loop, there 


are re-radiated fields in addition to the incident field. The total field, 


therefore, is 
B- Bs. Bs (3. 1) 


where B is the total field, B the incident field, and B’ the re-ra- 


diated field. From Maxwell's curl equation, 


= 0B 
VxE = a (3. 2) 


wt 
Assuming that the incoming signal has e! time dependence, the 


curl equation becomes 
VxE = -jwB = -jkcB (3. 3) 


where 


_ 27 _ w 
es Win 


Integrating both sides over loop, this becomes 


ff vxE. aS = -jk ff cB: dS (3, 4) 


3a an ad i hl ae 
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If the loop area is very small, then the magnetic field in 
the loop is uniform with magnitude By . The integral is 


e = -jke ff B- dS = -jkc By A (3.6) 


Where e is the induced emf; and where A is the area of a loop. De- 


fining an effective height of a loop antenna as 


ry 


h, = -jkA (3. 7) 


the low frequency input admittance of a loop with a constant current 
(Ref. 12) is 

jou ae ‘ 
Yo = [6z,ae +4" 66 di - | (3.8) 


An unloaded magnetic sensitivity K, is defined as 


K, = oM,/ (Ref. 10) (3.9) 
Now 
2 ae ae \ 
e = I)/ Yo = -jkA cBy (3. 10) 
Ih = (-jkA) Yo cBy = hy Yo cBy = AKp (cB) (3. 11) 
Ip , therefore, is a measure of the normal component of the magnetic 


field at the center of the loop, independently of the electric field. The 


loop sensitivity for the electric field, Kp , is a function of the size 


of the loop. 


58 


Let the current in the loop due to the incident fields consist 


of the zero and the first phase- sequence currents, i.e., 


I(s) = 1%(s) + {)(5) 


(3. 12) 


is +2/2) = 1%s)- 1%) 


Fig. 3.1 Rectangular current loop 


The first phase current 11)(,) depends directly upon the electric field 


(Rei, 10). Let the loop lie on the x-y plane and the center of the loop 
(1) 


at the origin, then I ‘(s) can be further broken up into 


4) 2 7), 1) (3. 15) 
x *"y 


corresponding to BY and eS . Therefore, the current at s = 0 
0 0 
is 


(0) = h YE (4) 


For a short dipole, hy = 419 where (¢ is the length of the antenna. 


Se 


59 


eae ae. ee Se ee eee 


Then, 
| (1) (i) 
r"(0) = AKL E 3. 15 
; y (0) > AKER (3. 15) 
where 
i 
Ky = Yhap / (3. 16) 
For a magnetic loop with a diameter D, the effective height hy becomes: 
h, = -ikm D*/4 (3. 17) 


and the magnetic sensitivity defined in Eq. 3.9 becomes 


sn ae 7D (3. 18) 


B D2 
Ang (@ - 3.52 + 13.0 2-) 
r»7/ 


where 


2 = 21n 
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2D 
a 
a = radius of wire loop 


The current induced on the portion of the loop parallel to 
the incident field is from the electric sensitivity. 


i D 


I(g) = : 
is) mye ja, 


cos ¢ (3.19) 
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where 7 is the characteristic impedance of the medium, and ay is 
an expansion coefficient given by Storer (Ref. 13). For D < 0.1, 


the electric sensitivity reduces to 


Kp = (j 2n° D*/2*)/ ny (Q- 3.52)(1-9.8D7/a#) (3.20) 


Fig. 3.2. Circular loop prove 


Also, it can be shown that for a loop witha diameter D<< 0.03A 


(3. 21) 


From the dependences of electric and magnetic sensitivities 
upon the loop size, and since the loop current is directly proportional 
to the sensitivities Kp and Kp , the loop diameter should be as small 
as possible to minimize the effect of olectric field upon the total current 


measured. The lower limit, however, is dictated by the fabrication 


technique and the receiver sensitivity. 


pee al et a_i a I 
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3.3 Experimental Procedures 

The current measurement on the antenna was obtained by 
detecting the magnetic field locally along the anten:.a with a circular 
probe that is sensitive to magnetic fields. Initia!ly, the experiment 
was carried out on a model at a frequency range between 5 MHz through 
30 MHz. The antenna used for the impedance measurement was 2.5 
meters long, which represents a quarter wavelength at 30 MHz. A 


one-quarter-inch copper tubing was used an an antenna, and the current 


Progrery teres 


along the antenna was detected with a traveling probe moved along an 
axial slit cut inthe antenna. This technique was first developed by 
ing et al for measurement of antenna current distributions. 
Since the object of this study was to examine the effects 
of the size and location of a finite-size ground plane, 2 systematic 
variation of ground plane diameters and their locations above the 
natural ground, up to a half wavelength at the highest frequency, was 
planned. In the frequency band chosen, this size represented a ground 
plane diameter and its location up to 5 meters. Due to the bulky physi- 
cal size of the full scale system, a smaller scale model wa« adopted. 
The scale models used were the same ones used in the 
radiation pattern studies. For current distribution measurements, 
the entire system was put into an anechoic chamber that is designed 
to operate above 100 MHz. The detecting probe for the local magnetic 


field along the antenna has a loop diameter much smaller than a wave- 


length, even at the nigher frequency of 200 MHz through 1200 MHz. 


Ae eT ae eee 
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The probe was suspended from the top of the anechoic chamber to the 
current element and its movement was synchronized using a Scientific- 
Atlanta recorder. The detected signal was then fed into a recorder and 
the result recorded. 

The probe, the anechoic chamber, and the equipment used 
for this measurement are shown in Fig. 3.3 through 3.5. The block 
diagrai: showing the current measnrement set-up is shown in Fig. 3. 6. 

Detection of a current along a current element with an 
external probe, as shown in our experiment, was accepted fully by 
many people. Since the near field may be disturbed by the presence 
of the probe, and its lead, a probe traveling inside an antenna with 
only that part of the probe which is used for detecting the local mag- 
netic field exposed through a slit was the optimum choice. However, 
if the probe is very small, compared with a wavelength, and if the 
lead from the probe to the recorder is well shielded from any stray 


pick-up, the effect of the external probe upon the antenna current was 


minimal. The lead from the probe to the recorder was a solid copper 
outer-conductor coaxial line (Micro Coax cabie) with an outer diameter 
of 0.5 mm. The probe lead was oriented perpendicular to the antenna 
current element, to minimize the disturbance of the near-zone electric 
field parallel to the current element. 
To test the accuracy of this type of current measurement, 
a quarter- wavelength monopole was mounted over a ground plane, that . 


simulated natural, infinite ground. The ground plane was a 5 wavelength 
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Fig. 3.3. Current probe 


SS ee Re ee ee ee 


s re 


64 


Fig. 3.4. Current measurement set up 


65 


Fig. 3.5. Anechoic chamber 
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Scientific- Atlanta 
recorder 


Probe 
movement 
control 


GR 120a 
Unit 
Oscillator 


Absorber 


GR 1203 
Unit 

Power 

Supply 


Fig. 3.6. Block diagram for current measurement 
set-up 
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67 
in diameter aluminum ground plane. The test frequency was 1200 MHz, 
where 6. 25 cm is a quarter-wavelength. The measured current dis- 
tribution at this frequency was plotted against an ©xact sinusoidal dis- 
tribution in Fig. 3.7. This measurement agrees very closely with the 
King measurement of a current on a cylindrical antenna with a probe 
traveling inside the cylindrical antenna with a slit cut along the axial 


direction. 


3.4 Measurement Results 


ad eatatieibitniteh Witibed “alsin 


The variables for the current measurements were the same 
as those for the radiation pattern measurement. In Fig. 3.8, a plot 
of the theoretical current distribution of a monopole over an infinitely 
large conducting ground is made for 6 different frequencies. “The an- 


tenna was a quarter-wavelength at the upper end cf the frequency, 


a al cl i dl ll 


namely 30 MHz, and as frequency was reduced to 5 MHz in 4 steps, 
the resulting antenna current distribution was rlotted in order to 
compare these currents with those on the same length of antenna 
over a finite ground plane located various distances away from the 

large ground. In Figs. 3.9 through 3.14 at each frequency the 

diameter of a small ground plane was fixed at 2 given fraction of 

wavelength and locations of this ground plane were varied in four 

steps. Figures 3.15 through 3.20 show the effect of ground 


plane diameter when located at a fixed distance away from the large 


ground. 
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Fig. 3.7. Comparison of a sine curve with a current distribution of 
a monopole over a large ground plane 
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Fig. 3. 8. 
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h = 2.5 Meters 


Infinite Cround Plane 


h = "4 at 30 MHz 
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Normalized Antenna Current 


T'eoretical current distribution of a monopole 
on an infinite grove ‘ane, 
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The same measurement was performed with a hemispherical 
ground plane, in order to study te validity of a sinusoidal current as- 
sumed on the artenna for a theoretical study which was carried out in 
Chapter 5. These measurement results are plotted in Figs. 3.21 
through 3. 22. 

Only amplitudes of the current on the antenna have been 
measured, and the maximum amplitude for each set of measurements 
was normalized to 1. 

The theory developed in Chapter 5 for an antenna on a hemi- 
spherical ground plane began with an assumption that the current dis- 
tribution on the antenna is in the sinusoidal ferm. It was further as- 
sumed that the antenna current was independent of the induced current 
on the ground plane surface. Accuracy of the theoretical results largely 
depend upon the accuracy of this assumption. From the measurement 
results shown in Fig. 3.21 through 3. 22, it can be seen that at frequen- 
cies where the antenna is near a quarter-wavelength long, the assump- 
tion of a sinusoidal current distribution was fairly accurate. However, 
at frequencies below 20 MHz, the currents on the antenna were modi- 
fied appreciably due to interaction between the antenna current and the 
ground plane current. At 5 MHz, in particular, the current distri- 
bution on the antenna was almost constant in magnitude, similar to the 
current distribution normally assumed for a Hertzian dipole. Origi- 


nal assumption of a sinusoidal current disty ibution should, on the other 


hand, provide a form ciose to a triangular distribution at 5 MHz 
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for an antenna resonating at 30 MHz. This type of distribution is simi- 
lar to the current distribution of so-called Abraham dipole. Assuming 
that the physical length of these two antennas were the same, their 
effective heights differ by a factor of 2. Therefore, a constant current 
distribution gives a radiation resistance as large as 4 times that of a 
triangular current distribution provided that the amplitude of the tri- 
angular currents are the same. 

A monopole with a finite disc antenna shows a better agree- 
ment between the assumed Sinusoidal current distribution and the ex- 
perimentally measured distribution at both the high and low ends of 
frequency band. Howevcr, at the midband, the current distributions 
show lower amplitudes compared with a sinusoidal distribution. Also 
noticed was that for a given size of disc diameters, the changes in the 
location of the ground plane height (a) affects the current distribution 
more than the changing ground plane diameters at a fixed position. It 
was also noticed that when the ground plane diameter was less than 
1.25 meters, the current distribution was such that the amplitudes were 


bigger than those of sinusoidal distribution. 
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Fig.3.18. Current distribution ona monopole antenna with a 
ground plane of various diameters at a given loca- 
tion with respect to an infinite ground at 15 MHz 
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Fig. 3.20. Current distribution on a monopole antenna with a 
ground plane of various diameters at a given loca- 
tion with respect to an infinite ground at 7.5 MHz 
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a hemispherical ground plane of various sizes 
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CHAPTER IV 


RADIATION PATTERN MEASUREMENTS 


4.1 Introduction 

The important parameters of any antenna are its radiation 
patterns and its input impedance. In this chapter, the experimental 
results of the radiation pattern measurements are given. The arrange- 
ment of the experimental results is intended to show the effect of the 
ground plane size and its location with respect to the infinitely large 
natural ground upon the radiation properties of the antenna system. 
Certain approximations and assumptions had to be made because of the 
practical limitations that confront this type of experiment, and they are 


discussed. 


4.2 Measurement Proble ms 

The input impedance measurements discussed in Chapter II 
were pe~formed with a monopole a quarter wavelength long at 30 MHz 
and & rand plane whose diameter varied from 0.625 meter up to 5 
meters. The experiment was conducted in the frequency band of 5 MHz 
through 30 MHz. In performing the radia‘ion pattern measurements, 
several assumptions and approximations were made. The antenna test 
range was originally designed for high frequency operation, mainly be- 
cause the required physical size and distance between the transmitting 


and receiving antenna were more practically realizable in the higher 
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frequency ranges. The linear dimensions of the original antenna and 


ground plane were reduced by a factor of 40 so that the operating test 
frequency was between 200 MHz to 1200 MHz. The electromagnetic 
scaling and conditions under which this is valid have been discussed in 
detail in Chapter 2. 

In the original antenna and ground plane measurement, the 
infinite ground below a finite ground plane was the natural ground with 
a certain electrical conductivity (0 ). One of the conditions for scaling 
the electromagnetic system is increasing the conductivity of the scale 
model by a scale factor. In this case, the large ground plane simulat- 
ing the infinite natural ground must have a conductivity of 40 0 . Also, 
the antenna and the ground plane for the scaled model should have been 
fabricated with a material whose conductivity was 40 times higher than 
the copper and aluminum used for the monopole and ground plane, 
respectively. 

For practical reasons, copper was used for both the antenna 
and for the finite ground plane. Copper mesh screen was used for simu- 
lating the infinite natural ground of the reduced model. The copper 
mesh screen used was in the form of a rectangular screen 5 feet by 
5 feet in size. This represents 6 wavelengths on one side at the upper 
end of the frequency band and 1 wavelength long for each side at the 
lower end of the band. 

Due to the finite size of this simulated infinite-ground- 


plane, the radiation patterns, in general, tend to have their maximum 
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radiation intensity off the hovizontal axis, where ordinary polar angle 
(8) in the spherical coordinate system is 90 degrees. 

A quarter wavelength monopole on an infinite conducting 
ground plane has a fadiation pattern in the form of sin@, which gives 
the maximum intensity of the radiation field in the direction of @ = 90°. 
Leitner and Spence (Ref. 3) have showed,theoretically, the dependence 
of this pattern on ground plane size. Fora circuiar ground plane of 
radius a, the far-field radiation patterns were plotted for ka = 3, 4, 5 
and 6.5. They all have maximum intensity off the horizontal axis 
where @ = 90°, 

Therefore, the results obtained from this experiment are 
expected to differ somewhat from the theoretical results where the 


ground plane was assumed to be infinite 


4.3 Experimental Technique 

The radiation patterns obtained in this experiment are taken 
in an x-z plane as a function of polar angle @, as shown in Fig. 4.1. 
When the antenna and the ground plane are oriented in this way, the 
axis of the monopole coincides with the Z axis of a rectangular coordi- 
nate system, the antenna is vertically polarized. This pattern is some- 
times called the "E plane" pattern where Eq: with @ constant asa 


function of @, is the quantity actually measured. 


The antenna range must satisfy a number of important con- 


ditions before it can be used for antenna radiation pattern measurements. 
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Fig. 4.1. Radiation patterns in the vertical plane 


The distance between the transmitting antenna and the test antenna must 
be sufficiently large to assure an accurate far-field radiation pattern. 
For an accurate far-field measurement, the antenna under test must 

be illuminated with a plane wave front. Since perfect plane wave is 


only possible at an infinite distance, some practical limits must be 


established as a guide to an acceptable distance between the transmitting 


and receiving test antenna. A common criteria is that the phase differ- 
ence between the center and the edge of the antenna under test be no 
greater than A/16 (Ref. 11), Then the minimum acceptable distance 


between the antennas is given by the equation 


R > 4a* (4.1) 
ae | 
where 
6<<R and 6<<a. (4. 2) 
The distance requirement is shown in Fig. 4.2. 
R? +a” = (R+6)? (4. 3) 
6 = (4. 4) 


The range used for this experiment was constructed with 
R= 20 meters. Since the radius of the finite ground plane, a, is as 
large as 4/4 


4a? 


es (4. 5) 


IA 
ha] >” 


where A varies from 1.5 meters to 0.25 meter from 200 MHz to 1200 
MHz, the distance R is more than enough to satisfy the condition (4. 1). 


Another possible source of error in the radiation pattern 


measurement is due to ground reflection, when both the transmitting 
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Transmitting Receiving 
Antenna Test Antenna 


Fig. 4.2. Phase difference between center and edge 
of the test antenna 
and receiving antennas are located close to the ground level. One way 
to avoid this problem is to utilize a narrow beam transmitting antenna. 
Another possibility is, of course, locating both antennas far above the 
ground level. 

The transmitting and receiving antennas for this experiment 
were located approximately 10 meters above the roof level of a three- 
story building. The transmitting antenna used was a wide-band, fre- 
quency independent, log-periodic sheet triangular-tooth structure 
designed to operate above 150 MHz. It was linearly polarized and had 


a beam width of approximately 60°. 


The test setup used for the experiment is shown in Fig. 4.3. 
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Fig. 4.3 Radiation pattern measurement Set-up, block diagram 


4.4 Measurement Results 


The >per ment performed here was designed to observe the 
effect of a finite yround plane upon the radiation cha:acteristics of a 
monopole displaced a certain distance above an infin:'* conducting 
ground. The gecmetrical arrangement of each component of this ° i- 
tenna system is shown in "Wig. 4. 4. 

The values used for a, h, and Z in the exper.ments were as 
follows. The physical length of h for Fig. 4. 4 was fixed so that it is 
one quarter wavelength at the upper end of the test frequencies, which 
ranged from 200 MHz threugh 1200 MHz for a flat disc ground plane. 
The spherical ground plane shown in Fig. 4. 4(b) also had a fixed length 
of antenna which was one quarter wavelength at the upper end of the fre- 
quency band ranging from 246 MHz through 1475 MHz. 

The frequencies chosen for the spherical ground plane were 
based mainly on the available size of copper half spheres. A quarter 
wavelength at 1475 MHz is 2 inches or 5.08 cm. A set of half-spheres 
available was 2 inches, 3 inches, 4 inches and 5 inches in diameter, 
and at the frequency of 1475 MHz, they correspond to 4/4, 3/8, »/2, 
and 54/8 in diameter. The radiation patterns are later studied asa 
function of ground plane size in terms of wavelength 2d. 

The experimental results are arranged first in Fig. 4.6 to 
see the effect of a ground plane size when it is located at a given height 
(a) above the infinite ground plane. The ground plane diameters were 


changed from 0.0156 meter to 0. 125 meter which simulates the 


Monopoie 
Antenna 


Circular Ground _, 


Plane a4 


Infinite Ground 


Monopole 
Antenna 


Semi-Spherical 
Ground Plane 


Infin] 
Infinite Ground 


(b) 


Fig. 4.4 Geometrical arrangement of an antenna and ground plane 


diameters of 0.625 meier to 5 meters of the original antenna system. 


These diameters correspond to 4/16 through A/2 at 1200 MHz and 
30 MHz. The changes in the diameter were made in four steps, with 
the size of the ground plane doubling its diameter each time. Thus, the 
diameters of the flat disc used for the experiment were /16, /8, 
4/4 and A/2. 

Second, in Fig. 4.10, the location (a) of the ground plane 
was changed for the given size of a ground plane diameter. The stens 
taken for different a were the same as discussed above for the variable 
size of ground plane, namely, a= A/16, 4/8, 4/4, and d/2. 

In Fig. 4.5, a monopole above a semi-spherical ground 
plane of various size has been tested for its radiation patterns. When 
the ground plane size becomes larger with respect to wavelengths, side 
lobes beceme roticeable. When ka=2, the side lobe level shows over 
60 percent xf the ..ajor lobe level. However, for smaller values of ka, 
side lobe effect is ue, ligible. Typical experimentally observed beam 
widths measured between 3db points range from approximately 30 degrees 
to 70 degrees. These yatterus are later compared with the theoretical 
results obtained in Chapter 5. 

Figures 4. §through 4. 13 show that when displacement be- 
tween the disc ground piane and a simulated infinite ground plane is less 


thar 4/2, major lohes are directed considerably off from the horizontal 
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position with a small side lobe along horizontal direction. This type 
of phenomena can be also observed for a center fed dipole of 1.5A in 
length (Ref. 14). This indicates that by using a monopole with a finite 
ground plane displaced a certain distance away from an infinite ground 
plane, the system behaves like a longer dipole antenna. A finite ground 
plane and a gap between the finite and infinite ground plane simulates an 
extension of an antenna length. It is also noticed that the radiation pat- 
terns generally are more sensitive to the gap distance a than the ground 
plane diameter D. In particular, when the gap distance a = A/4, 
the pattern shows a considerable difference fr »m other patterns. 

In summing up the results of radiation pattern studies, it 
can be concluded that the finite ground plane used in conjunction with 
an infinite ground below causes radiation patterns mainly to be directed 
off from the horizontal direction due to an equivalent antenna arm created 
by the finite ground plane and the distance to the infinite ground. Although 
minor lobes are present in most cases, their relative levels, compared 
with major lobes, are insignificant. Also, it shows that t-\e distance 
a is more significant for the radiation patterns than the diameter of 
a finite ground plane D. Finally, there seems to be sc ne similarities 
in radiation patterns between the monopole with hemispherical ground 
plane and the sleeve dipole (Ref. 15). However, a detailed comparison 
is not possible because very few corresponding dimensions between two 


antennas are available. 


Fig. 4.5. Radiation patterns 1G, [? fora monopole with a hemispherical ground plane of 


various diameters 
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Fig. 4.6. Radiation patterns lee? for various ground plane 


diameters (D) with a = 0, 62"r 
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Fig. 4,7. Radiation patterns |E | for various ground plane 


diameters (D) with a = 1.25m 
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Fig. 4.8. Radiation patterns lE, |? for various zround plane 
diameters (D) with a = 2.5m. 
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Fig. 4.9. Radiation patterns IE, | * for various ground plane 
diameters (D) with a = 5.0m. 
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Fig. 4.10. Radiation patterns |E 9 1? for various ground 
plane locations (a) withD = 0, 625m 
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Fig.4.11. Radiation patterns IE 1° for various zround plane 


locations (a) with D = 1,25m 
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Fig. 4, 12,Radiation patterns |£ | ® for various ground plane 
locations (a) with D = 2,50m 
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CHAPTER V 
THEORETICAL ANALYSIS: 


5.1 Introduction 

In this chapter, the theoretical model analyzed is a 
monopote antenna on a hemispherical ground plane that is, in turn, 
above an infinitely large conducting ground. Far zone electromag- 
netic fields are computed as a function of hemisphere radius and an- 


tenna length. The model chosen for the analysis is shown in the Fig. 5. 1. 


Monopole Antenna 


ini Hemispher ical 
Infinite 
Conducting Ground Plane 
Plane 


(a) 
Monopole 
Antenna 


Image 
Anten 


(b) 
(c) 


Fig. 5.1. Theoretical model of an antenna system 
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Using image technique, the criginal model is divided into 
two parts as shown in the Fig. 5. 1(b) and(c). Each configuration is 
separately studied and later the far-zone fields due to both antennas 
and the induced currents on the spherical surfaces are superimposed 
to obtain the totai fields. 

The Section 5.2 deals with a development of a general far- 
zone fields expression from a given source configuration, and the Sec- 
tion 5. 3 includes a far-zone field calculation due to two monopoles 
separated by the diameter of the sphere. The Section 5. 4 computes 
the induced current on the spherical surface due to the monopoie and 
its image. And, the Section 5. 5 shows the far-zone field expressions 
due to the induced surface currents on the sphere, and the total far- 
zone fields due to the antennas and the surface currents on the hemi- 
sphere. Some of the results from the numerical calculation in the 
forms of radiation patterns are given. Section 5. 6 contains results of 
radiation resistance theoretically obtained for a monopole with a hemi- 


spherical ground plane. 


0.2 General Far-Zone Field Expressions 


0.2.1 Classical Formulation. Let us posiulate a new 


form of Maxwell's equations for harmonically oscillating fields: 
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4 
* 
2) 
N 
ee 
€ 
i= 

o 
x 

( 
7 
4 
al 
tt 


“JOP, 
(5. 1) 
VxH = J +jwe E V-J = -jwp 


where 


J Ps is a fictitious "magnetic current density" 


and 
Pm is a fictitious "magnetic charge density" 
From the above equations, the following equations are 
derived. 
V-(€oE) = p V+ (uoH) = 2, (5.2) 


Now define the electric and magnetic fields as sums of two parts, 
i.€., 
E = E'+E" 


2 (5. 3) 
Hi 


u 
= 
+ 
et 


where E' and H' correspond to the fields arising from the actual 


current J when De = 0 and thus satisfying the equations 


VXE' = ~ jw oH V+ (€ 9’) = p 


(5. 4) 


VxH' = J + jwe)k' Vv: 


—~ 
_o 
(on) 
mi 
" 
So 
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Similarly, E” and H" satisfy the equations 


Vx" = -jou,H" - Jin v= (€9E") 


VxH" = jwe Ee" V - (uoH") 


"1 


(5. 5) 


For the set of equations 5.4 , a vector potential A 


and a scalar potential @ , can be defined as 
HoH = VxA 


and 


where 


E'+jwA = -V® 


Then the following two differential equations are obtained. 


VA+KA = = Upd 


Vor+rko = -& 
*0 


where 


(5. 6) 


(5.7) 


(5. 8) 


(5.9) 
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The solutions of these differential equations are 


A= uy If GCS ) 50") av’ 
y' 
7 (5. 10) 
1 R pt t 
@== fffGi—) plR')av 
Oy 
where 
Ry eikIR - R'| 
OS) = ar iR-R (5. 11) 
and 


R and R'are the radial distance from the origin of the 
coordinate system to the observation points and the source 


points, respectively. 


Fields E' and #' can be written as 


E' = -jw[A+1 viv A)] (5. 12) 
k 
and 
H' = i VxA 
Ho (5. 13) 


Similarly, for the set of equations 5.5 , a new vector 


potential An and a scalar potential ¢ can be defined 


Ca gue Bade et eae ch ee Eros yt eis eee nS POE ESE I EE ede ER Oe Reece | emg ee mee am oa Me ar acre, ogee en eRe eee gh hay Ree eA SBME EMS 
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€.E" = -VxA 
0 m 
-jw = V'‘A 
Jon 9¢ 9 o” m m 


The resultant differential equations are 
ze 27 __ Ft 
Vv? At AL = ~€od 


and 


The solutions of these equations are 


>! 


we sHN0(B) 


spelt RY. (aay 
haa {J Jo() o,,vav 


The fields E" and H" can then be expressed 


2 1 as 
BE" = -— VxA 
a) m 


mute . os A b.. ck 
H" = -jw| Ant V(V A)! 


I ae ee WE ae ee, ee eS RT 


(5. 14) 


(5. 15) 


(5. 16) 


(5. 17) 


(5. 18) 


(5. 19) 
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Therefore, the total E and H fields due to both electric 


and magnetic sources can be written in terms of vector potentials as 


ae § ae. ess es | x 1 < 
& 1 moo os . ete VV. 
E = E'+5" 2 Eee V(VvV - A)] eo x A (5.20) 
r a. l Ts 1 
H = H' H" = — V - —V ° 
ig er erm ml AGRI) 
In the region where J = J = 0, these equations are 
simplified to 
as _ jw = l = 
E = Ne VANE A ONS AS (5. 22) 
0 
He  vxvxd_ ++ xa (5. 23) 


ke? eG) 
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Y 
A 


P(R, 8, 9) 
Observation point 


P'(R’, 8; 9) 
Source 
point 


Radiator 


Fig. 5.2. Coordinate system used to derive far-zone electromagnetic 
fields 
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Consider the coordinate system shown .above in the Fig. 5. 2 for far- 


zone field R >> R' and IR-R'| = R- R'cos a = IRI - (R'- 


“A “A “A 
If the direction of R' is denoted by @, @ where R, and R' are 


the unit vectors, then 
cos a = sin@ sin @'cos (¢'- ¢) + cos @ cos @' 


With the approximation shown above 


IR- R'l = R- R' cos a 
and 
‘i e IkR 
a1 0 ~/51,_jJKR' cos @,,,, 
A> pe SS STR dV 
Similarly 
€ @ Kh kR 
- <0 = 1 JkR'cos a ,,,, 
7 Sa JSJF (Re dV 
Then 
ye ier ; 
7 EM e JKR 7 
Pane edge MO) 


(5. 24) 


(5. 25) 


(5. 26) 


(5.27) 


(5. 28) 
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where 


Ho, 9) = {fT ekR’ - Bay. 


(5. 29) 


a 5D URY ans 
Mio, 9) = ff {—2— elke R ay: 
9 
N andM are only functions of 6 and ¢. 
N and M are called the electric and magnetic radiation vector, respect- 
ively. 
In general, 


N = NAR + NAO + NB rr Nak + Ny, 


where N is a transverse component. 


* 'R sino 


oN n 
1 20 oes Bee, 5. 31) 
Th (sin @ Ng) R ( 
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in spherical coordinate system, if R> > A, then all terms except 


jkN, and jkN 


can be. neglected. Therefore, 


¢ 
ager A Lge 
vx = 7 jk [1,0 - Ni] = -jk-—“S 
Since 
= ue eR = - 
VxA = Vx — iE N(6, ®) = VxON 
= -jk®(R x N) 
where 
i e IkR 
o = 0 
47R 
VxVxA = Vx(VxA) 


where we cantake Vx A aS a new vector (W) equivalent to N above. 


VxVxA 


VxW = (-jk) @x W) 


ke? (“s 


ie jkR 


4nR 


aN /\ 2 Ky 
(Rx(RXN)) =k — 


where 


~jkR 
are -N ui®xa e 
E = jkn,[ N. +R x M, | ih 
secs (5. 34) 
a _ ostta 2 ed lA Cs 
H = jk [M,+RxN, GR nee 


These expressions show that only the transvers? com,,7nents 
of Mand N enter into E and H far-zone expressions. 

5. 2.2 Stratton-Chu Integral Formulation. In the previous 
section, a general far-field expression of E and H have been developed 
from Maxwell's equations with a given current or source distribution. 

In this section, the same far-field expressions of electromagnetic fields 
are developed, using Stratton-Chu integral formulations. 

Let P and Q be two vector functions of position with 

: the proper continuity, then 
ee (Q:VxVxP-P-VxVxQ) av 
(5. 35) 
= fT REV RA OUR E YES 


where S is a regular surface bounding the volume V. The above 


integral equation is the generalized Stoke's theorem. 


ee eet ner ge me mmo Fee 
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If the field vectors E and H are assumed to be e it 


time dependent, the Maxwell's equations can be written as 


VxE- jwuH = -J 


The medium inside volume V is considered to be homo- 
geneous and isotropic. The quantities a and Py are the fictitious 


densities of magnetic current and magnetic charge. 


Fig. 5.3 Notations for Stoke's theorem 
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We consider a volume V_ bounded by the closed surface 


Sy Sy) Sy ‘aie sy as shown in Fig. 5.3. The notation S in Eq. 5.35 
means a sum of Sy Ss are s: 


From Maxwell's equation, expressed in Eqs. 5.36, the pair 


of vector Helmholtz equations can be derived 


VxVxXE-KE = jopd - Vxd 


(5. 37) 
7xVxH-k’H = joes + VxJd 
Where currents and charges of each are related by the continuity 
equations 
3 
V-id- jwp = 0 
(5. 38) 
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and 
ke = wpe 
In Eq. 5.35, let P = # and Q = ¢@ where @ is an 
; arbitrary unit vector and ¢ = elkt / r. Distance r is measured from 


the element at (x', y', z') to the point of observation (x, y, 2); i.e., 


r = (x= x’? tly-y)? +(z2-z'? = IR- Rl 


(5. 39) 


a 1 
R = (x? +y* +27)? and R' = (x” +y” +2")? 


Then, 
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= AN 
VxQ = Véxa 


(5. 40) 


at las Nn 
VxVxQ = ak’¢+V(a- V9) 


Vr VxP = E+ jopJ-Vxd 


Substituting Eqs. 5. 40 into 5.35 gives 
= - = A 
SSS [va (kK E+ jwpI- VxJd_)-E- (akg+v(a v9) ] av 


= ff [Ex Véxa-SaxVxE]-dS (5. 41) 


Since 
_ Fan AN oo SN a 
E-V(a° V¢) = V-e(a* VAE- (a° VO)V-E (5. 42) 


and 


= 1 
V°-E =p 


[Sf v-@-veEav = ff (@-vo)E- a8 
V S 


=@. ff (mE) vedas (5, 43) 
s 


Equation 5, 41 becomes 


SIS [iwud ¢ - VxJ_o + t pVv¢] dav 


= ff [jw p(x AN + (Nx E) x VG +(- E) VO - Axd mi?! 98 
Ss 
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where 


Notice that 


JS{(@xvdxa)- fds = ff ix (Ex v¢)-aas 


and 


Sf (@axVxE)- nas = Jf (ax (jon -T_) - n aS 


it 


- ff jou(x fs. aas- ff ixd 6+ ads 
Equation 5. 41 can be further reduced to 
SSS GouIe-F_xve+4 p V¢) av 
V m € 


= ff 'jwu(fix B) + (ixE) x Vd + (fi. E) V6] ds 
s 
(5. 44) 


with the help of identities: 
SSfV x G9) av = Sffvoxd, ave Sffovxd, av 


and (5, 45) 


SJf ¥x@,p)av = ff ixd gas 
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Thus 


~ff[JovxT ev = -ffaxT, eas- ff fT xveav (5.46) 


In Fig. 5.3, the surface S includes S1, So» ait s, and 


a volume surface So surrounding the observation point (x, y, 2). 


Since 


- 
S—) = (-- iw) -F 


Vera Ve 0 (5.47) 


where To is a radial unit vector pointed toward the center 
of the spherical surface S,. Thus, fy = for the surface S, . 
Evaluating the surface integral of the Eq. 5. 44 over the surface 


Sy » notice that 


(nx E)x@+ (@ E)n = -Px @xE) + (@ E)n 


“ME 1) + Enh) + (@- ED (5. 48) 
= £ 


Also, when the radius of the sphere Ty 


is reduced to zero, the only terms remaining from the surface integral 


represented by So 


are terms involving 1/ry* 
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fJ [ @x ) x Vo + juu(ax Ho + (n= E) vo | as 


S 
2 a jkr ee! as jkr 
=f f (@xE)xt oN + By Bye, £2 + of 1 ) 
0 0 0 : r : r” r” 
0 0 0 
ry" sin@ do do (5. 49) 
where 0/ ) represents other terms involving A 
*9 To 


As r, approaches zero, the integral tecories 


rc "IN - IN cm IN an 
JJ U(nx E) x Vo + jwu(nx H)o+(n- E)Ve)as' 


= 4nE(x, y,z) (5. 50) 
Equation 5.50 then becomes 


E(x,y,2) = Ph eRe dak Vo +— p Ve) av" 


- JJ [jwu(Tx H)o x (nx EB) x Vo + (4. B) Vo] as’ 


S,, 85, ..-8) (5,51) 


If the volume does not contain any sources, Eq. 5.51 further 


reduces to 
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E(x, y, z) ee 7 Q 
1°~2?°° "nn 


[jw ux H)¢ + (ix E) x Vg +(t- E) Vg] ds’ (5. 52) 


Also, if all currents and charges can be enclosed within a sphere of 
finite radius, then the field is regular at infinity and either side of S 
may be chosen as its interior. The surface integral of Eq. 5.51 rep- 
rescnts the contribution of sources located outside S. If, therefore, 
S recedes to infinity, the contribution from those sources vanishes. 


Discarding the fictitious magnetic charges Jn , Eq. 5.51 becomes 


Bix, vs 2) = gz SSS [wu to+—2p v9] av (5. 53) 


Equation 5. 52 is the field that would be produced from surface S by 
electric current density K, by magnetic current density ae , and by 
surface electric charge of density Oe 5 where 
—_< oe oo ana be 
@xH =-K, ®xE = Ky % = et E (5. 54) 


If Eqs. 5.52 and 5. 53 are applied to the far-zone where 


jkr 
g = el /r, then 


QikR - iKR - Re 
¢=——_———_ (5. 56) 
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and Eqs. 5.52 and 5.53 can be writt2n as 


E(x, y,z) = x ff [jwpe(fix He + (ixE)x V¢ J, dS! 
$1,59--- 8, (5. 56) 


and if 8), ) wr s. are perfectly conducting boundary, then ix E = 6 


Therefore, 


E(x, y,z) = = ff [jwu(n x H)¢ ], ds’ (5.57) 
4n S..s s t 
159, --- 8, 


where subscript t signifies transverse component to the radial direc- 
tion. The proof of these equations 5.56 and 5.57 is given in Appendix 
B. Of course, Eq. 5.56 is equivalent to Eq. 5.34 in Section 5. 2. 1. 


Similarly, it can also be shown that in the far-zone, 


A(x, y,2) = = ff [jwe(fix E)¢ + (fix H)x V9], ds" 
51) So. me 5, 
(5.58) 


5.3 Far-Zone Field Expressions for Two Linear Antennas 


Consider an antenna and its image with a sinusoidal current 


distribution as shown in Fig. 5. 4. 
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Fig. 5.4. Collinear dipole 


Let the currents on antennas A and B be 


Tle") — ‘ _ ot) 1 
I(z') = Tax 510 kid “)z foraca<d 


(5. 59) 


I, Sin k(d + 2") 2 for -d < 2! <-a 
The Antenna B is the image of the Antenna A. 


Let us find the expression for the electric radiation vector 


N. Since we are only considering the electric source 
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. Dt “ 
= JJ JT Ry ekR'Ry,. (5, 29) 
where 

as ran 
R'- R = R'cos a = R' cos 8 cos @' + Sin 6 cos (9 - 9’) 
™N 
R = Unit vector in radial direction 


In the problem considered only line sources are present. Therefore, 


= pei TERR 
N= fie! ae (5. 60) 
where R' = 2'2 and R'- R = 3! cos 6 for antenna A 
and R' =-2'2 and R’- R = +Z' cos 6 for antenna B 
d ike’ cos, A 
N= f T, Sin k (d- 2") elX@" Cos ‘ Ee sin k (d + z') elK2' cos 0, , 
a -d 
(5. 61) 
Changing the variable for the Second integral by z' = -2"", 
the integral becomes 
-a qd 
ff m a k(d+ z'tve ,}k2' cos 8 ae = f In Sin k(d - 2") e” Ik2" cos @ dz" 
= a 
(5. 62) 
Therefore, 


d ak ees 
N = Pa sin k (d- 2'){ elk? Swe are jkz" cos / dz' 


d 
= 2In | Sin k(d - z') cos (kz' cos 8) dz' 


¢ (5, 63) 
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Since 


; { sin (a+b) + sin {a-b)] = sinacosb 


d 
N = et | {sin puta 2°) + cos 6] + sin [ k (d- 2') - kz’ “ 
a 


d 
= 1, f sin [ kd + kz'(gas-1)] + sin [ kd - kz'(1 + cos @)] dz" 
(5. 64) 


Also, using 


f Sin (a - bx)dx = + cos(a - bx) 


1 wl. 1 
N, = In econ) cos[ kd - kz'(1 - cos @)] * Eis cos 0) 


d 
cos [kd - kz'(1 + cos 6)}] 
a 


es l l 
= 1 R= cos 0) cos [ kd cos @ } +e Wa cos 6) cos[ kd cos 0] 


cos [ kd - ka(l - cos @)] - k cos [ kd - ka(1 + cos af 


l l 
~ kl - cos 0) 1+ cos 0 
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N_ = 1_ cos (kd cos 9) fo + : J 
Zz m k(l- cos 0) ~ k(I+ cos 6) 


-y eos k(d - a) + ka cos @ , £08 k(d - a) - ka cos 0 
kil - cos k(@1 + cos 


2 cos (kd cos Q) I 
Sree eS, es fa + cos 9) cos[ k(d- a) + ka cos 9] 
k sin? @ k sin? @ 


+ (1- cos @) cos [{ k(d- a) - ka cos 0| 


Using trigonometric identities, Eq. 5.65 can be written as 


2 cos (kd cos Q) aI 
\y a ns { cos k(d - a) cos (ka cos 9) 
k sin? @ k sin? @ 
- cos @ sin k(d - a) sin (ka cos @)] 
2] 


k sin? @ 


| cos(kd cos @) - cos k(d - a) cos(ka cos @) + cos @ sin k(d - a) 


sin(ka cos @)] 
(5. 66) 


When a = 0, the expression reduces to the radiation vector of a dipole 


which is 


2I 
No ee | cos(kd cos @) - cos kd] 
k sin* 0 
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The @ component of N is 


Nog = -N, sin® = in { cos(kd cos @) - cos k(d - a) cos(ka cos @) 
+ cos @ sin k(d - a) sin(ka cos 9)] 
(5. 68) 
The far-zone electric field is, from Eq. 5. 34, 
-jkR 
E = jkn (-N) “aR (5. 34) 


when 


Therefore, the electric field due totwo monopoles, E q? becomes 


: -jkR 
E. = I In id EF: cos 9) - cos k(d - a) cos (ka cos Q) 


oq 2m sin @ 


, 0s @ sin k (d - a) Sin (ka cos 8) 


sin 9 


(5. 69) 


F(Q) = 
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The radiation patterns due to the two monopoles discussed 
above are the plots of angular dependent terms of E as a function of 


6. 


cos(kd cos 9) - cos k(d - a) cos(ka cos Q) + cos @ sin k(d - a) sin(ka cos 9) 


sin 0 (5.71) 
When 
a=0 
F(0)= cos(kd cos 9) - cos kd (5. 72) 
sin 0 


This expression is the angular term of a far-zone field expression of a 
dipole. 
5.4 Induced Current on a Spherical Surface 

5. 4. 1 Induced Current on a Conducting Sphere Excited by 
a Monopole. The induced current on a conducting sphere excited by a 
monopole antenna erected on its surface was treated by Papas and King 
(Ref. 16). The current along the antenna was assumed to have a form of 
sinusoidal distribution and was assumed to be independent of the current 
on the sphere. 

This solution is reviewed here, because the same method 


will be adopted to calculate the induced current on the spherical surface 


due to the image antenna located along the negative Z axis. 
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Xx 


Fig. 5.5. A monopole antenna above a spherical ground 
plane. 


From the general reciprocity theorem (Appendix A) 


ye eee 
ee a eae ee (5. 73) 

where Ey and Hy are the electromagnetic fields associated with a 

current density J, and E,, Hy with J, Also s_ is the closed surface 


of the volume of empty space bounded by the surface of the antenna and 
the sphere and an imaginary boundary at infinity. f is an external normal 
unit vector of this surface. 
The assumption made here is that the two sets of fields, 
Ey Hy and Ey, Hs are defined in the free space and are bounded 


by the same geometrical surfaces, but not necessarily satisfying the 


same physical properties. 
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Let fields E 4 and H, be caused by the actual currents 
in the antenna located at 6 = 0 anda < R < a+h as Shown in 
Fig. 5.4 and the current on the surface of the sphere. 


The boundary conditions for fields E, and H, are: 


n x E, = 0 on the spherical surface 5S. 


Ww wot 7) " 
t x Hy ~ “Ls 8, 


The current on antennas Ty, A and K.. on the sphere is 
assumed to be maintained by an appropriate generator or a distri- 
bution of generators in the antenna in such a way that no current 


exists inthe » direction of the antenna or on the sphere. 


Therefore, 


K, (8) = OX, (0) (5. 74) 


and 


2a 


J K,,(0) asin@dg = 1, (0) (5.75) 


Ss 
1, (9) is the total current crossing a parallel of latitude on 
the sphere. 
Also, the total current in the antenna at @ = 0 is in radial 


direction in the spherical coordinate system. Due to the rotational 
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symmetry about the axis 9 = 0, the line integral of the tangential 
component of H, around the surface of the antenna is equal to the 
tota! axial current 1, (8) in the antenna. 


Thus, 
- = 2u _ 
OMA eS Big tg ff 1, /d8 = Ty 


(cr the —) 


of the antenna (5. 76) 


The second set of fields E, and H, are not the actual 


icids that are related to the real currents iy in the antenna or on 


the sphere. They are, instead, fictitious fields that satisfy the same 
field equations as the first set and that must be defined over the 
same geometrical surfaces. Since fields E, and H, are not 
required to be defined on the boundaries with the same physical 
properties as set E, and Hy let the volume possessed by the 
antenna for Ey anc Hy be empty space for Ey and H,. Also, 

on the surface of the sphere, S., fields E, and Hy are required 
to possess the following prescribed form by some appropriate 


set of generators: 
29 = U sin ® P (cos 0) 


(5. 77) 
= 0 


The entire field in empty space, therefore, has the 


components Eop» Eo, and Hog which must satisfy the field 


equations and the prescribed boundary conditions. 
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Returning to Eq. 5. 73, 
JS fn @, xi, - Bx] dS = 0 (5.73) 


the surface S_ can be divided into two parts S m the surface of the 


sphere, and S,, the surface of the antenna. 


os ox = is es xo _ ces NN 
Sf (8 xB, - B x Hi] d& = Sta (xB) +£, - (n'x H,)} a8 
s 


+ ff (Hi, -@x 8) - 8, -@x #,); dS 
S 
A 
(5, 78) 
N F ran “~~ : yas ran 
Ow since n = <@ over Sa and since n = -R over § and 
s 


with the boundary conditions described above, each integral can be 


rewritten 
{5 at+h Qn 
{fH + @xE,) - i, - (nx E.) aS = f 
2 yo" 2 E,p dR { Hy, Ry d 
S, } Ria OR 2g “9 a0 
i ie 
~ E,, GR | H,,R, do 
R= 2R rs 1g 0 
(5. 79) 
The volume occupied by the antenna is empty space for 
H 


Therefore, for a thin antenna where the radius Ro goes to zero, 


29 ° 
the line integral 
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2a 
{ H, Ro os = 0 
Therefore, 
ath 
St H,- Ax E,)- i, (ix B)}as 2 =f Eon a(R) aR 
‘A 


(5. 80) 


The surface integral on the sphere becomes 


StH, (xB) +E,-@xi)] 6 = ff B,- K,, as 


(5, 81) 
because of the noundary condition nx E, = 0 on §,. 
{ 
From Egqa. 5..75 and 5. 77 
Kn f dQ 
{ff B+ Kj, 48 = J Egg 1,(@) a 
s 0 
s 
Therefore, 
Roe fot. 8 = a 
ffl 2: @ x Hy ~ EB x H))] dS = 0 (5. 73) 
s 
which implies that 
j ie 
E,, 1,.(0) ad® = Eoy 1, (R) dR 
6 20 ‘Is A 2R 1A (5, 82) 


This is the integral equation to determine 11,(9) knowing 


1 AfR) and the prescribed field Eo, will assume that the surface 


io; fe it gee at ee eee ee ake, ne ee ee ee ee er eee Be SO eT ee er ees ee cyial net SoS aaa eG ea ea 
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current 1, (9) takes the form of 


© 
1,(0) = 2 BP (cosé) (5. 83) 
n=0 


where the coefficients BL are to be evaluated from the integral 
equation. 

For the second set of fields E, and H,, the volume 
surrounded by the closed surface 83 Ss N and the boundary at in‘inity is 


empty. Furthermore, Eo Eo, and H,,, are the only non-zero 


R’ 2¢ 
components due to currents on the sphere. For TM fields of this type, 


the fields are found from the scalar wave equation 
2,U 2,uU 
Vv (R)+k (F) = 0 (5. 84) 


where 


The solution of the equation is 


(o 0) 
u = L AP (CoS 6) p (KR) (5. 85) 


where 


P (cos 6) is the Legendre polynomial of order m 
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(2) 
p(k) = fy 1 OR) = (eR) n) aR) 


Hr) : Hankel function of the 2nd kind 
(2) 
h m ER) : Spherical Hankel function of the 2nd kind 


The components of the fields are obtained from solution u as 


2 


E,, = [k? + ]u 
a aR? 
5 Lee 
Fog = 8 ‘GR oO Y (5. 86) 


~ wi (K)/_OU 
From Eq. 5. 85 and 5.86 , on the spherical surface R = a, 


E = ~&) Y A_P Nicos 6) "(ka) 
20 °»~C lf a uo mm m (5. 87) 


Also, we demand that Eo take on the spherical surface: 


Pi , cos Q) - Pea jXcos 6) 


En, = U sin @ P (cos 8) = U mal 


20 


(5. 88) 
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Equating the two equations 5.87 and 5. 88 leads to 


A oe a U 
n+l n+)" kp, , (a) 
(5. 89) 
¥ a U 
An-1 = 2n+l): k- p ka 


n-1l 


Therefore, 


es te E es aP,_ j(cos @) p, _ {(kR) ; aP, (C8 9) py , (KR) 7 
2R aR? (2n +i): k-p__ ,'(ka) (2n + 1) ~k- Py’ (ka) ; 1 (ka) 


f 
n- 1 


(5. 90) 
Assuming that the ant2nna current takes the sinusoidal form 


I, , (R) = Te sin k(d - R) (5. 91) 


it can be shown that (appendix c) 


a+h Ual p_ - l(kd) p (ka) p, (kd) 
f Beowd (Rida Sa a - cos kh n=l Ail 
s 2R 1A an+1 p, - 1' (ka) p.- 11 (ka) Ba os ka) 
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Also 


T Tt oO 
f Bog 1,,(0) ad@ J Eog 2), B,P, (cos 6) ado 


n=0 


(5. 93) 


Finally, the integral equation produces the result 


Tax Py. (kd) Pn. (ka) Pus (kd) 


BOo=—s |] - ‘s 
n 2 Pp, y (ka) cath p,- I' (ka) P, 1 (ka) 


The surface current I, gt) is now determined in terms of the 


coefficients B 


[o8) 
1,,(6) = d B LP (cos 6) (5. 83) 


5.4.2 Induced Current Due to the Image Antenna. In the 
Section 5. 4.1, the induced current on the spherical surface due to the 
sinusoidally distributed current on the antenna was found in the form of a 
modal current. The coefficients B of the infinite series that gives 
the surface current were evaluated in terms of known parameters. In 


this section, the contribution from the image antenna to the total surtace 


current on the spherical ground plane is obtained following the method 
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used in the previous section. 


Fig.5. 6. Image antenna with a Spherical] ground 
plane. 


The geometrical arrangement of the image antenna and the spherical 
ground plane is shown in Fig. 5.6. The center of the spnere is 


located at the origin of the Spherical coordinate system. 


SIRE a OF Ly PT Tae OF Ee a 
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The electromagnetic fields and the currents are all expressed 
in the same way as they were in the previous section except for the prime 
notation to distinguish the image antenna problem from the real one. 
From the reciprocity theorem, 

SY (E, xi, - E,x Hi): dd = 0 (5. 95) 
where surface S is now the closed surface enclosing the volume bounded 
by the spherical surface, the surface of the image antenna and a boundary 
at infinity. The surface normal unit vector ‘Y is the same as was in 
Section 5. 4. 1 on spherical surface S, . However, on the surface of 
the image antenna S A the unit vector is in the positive @ direction. 


The integral equation 5. 82 is now modified to become 


| re ee ath 
J Eo, 1,,(0) ade = - f Eon 1, (R) aR (5. 96) 


Let the surface current Is (@) onthe sphere due to the 
image antenna also take a modal form with coefficients of the infinite 


t 
series Bo: Then, 


© ' 
1,,(6) =), BLP, (cos @) (5. 97) 
n=0 


Following the steps shown in Section 5. 4. 1 


ae. aP,_ (cos Q)p (KR) (cos 0), See 


Piel 


Se CHEE Kp. _ yp tka) ~ set ace LS ta (5.98) 


n+l 
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Since 


-l 
Pp _ (cos Q@) = (-1)" atQ=7 


(5. 99) 
along the image antenna 
P+ 1'68 @) = (-h*! at @ =n 
Eog = U[k? +2 pre! etal yt Be 
n+1) Kp _ )' (ka) (2n+1)-k Py (ka) 


= 1} ope + ZY P _ \(kR) 7 52208) 
ak? (2n+1 kp 7, ka aa OT ic P 
ath Ual (ka) 
c t ' = a z 7 n ee Pn “] 
7 "aR nae = J cy" For a’ OR = (-) 2n +1 a lkay 
| Pui (ka) 
i - cos k(d - a) td net + cos k(d - a) P+) 
Paka) py, 1’ (ka) te oe 
Since 
ath 


Por T, (R)dR 3 oe Ey9'1,,'(9) adQ 
me (5. 102) 


OT SE PET REE ae OE TT NPE eee ey 


Pe Bat Sah tas, cial A as Bi 


2B, 
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where 
' : 
1 a*(R) = Le? Sin k(d - R) 
and 
00 
1,9) = ) B,' P.(cos 8) (5. 97) 
n=0 
Also 


P ‘(cos 0) - P "(cos @) 
Ey9' = Usin@ P (cosa) = u{ At} ___n-l 


= 5. 103 
2n+1 | ( ) 
Thus, 

I Boo! Y BP (cos@)-ad0 = B' U2a 

5 20 Zp non n (2n41) (5.104) 


Using the orthogonality relationship for the Legendre function 


em 
ae] 
3 


(u) Paw) dn = Ear 18m 


Pp, - a) p 


; , 2 pd Fess WPA a > 
= f Eo, 1 ,'(R) dR = (1) Tinax (p — (ka) cos (4 - a) >——a75) ea 


Pual k(d Pn (a) 
“By, tea + COS a5 .*(ka) 
Or 
i+ 
ee Cle ae Pa: 1) coc en Paes) Poy 1 
B= 2 Py. 1 (ka) Ph - "(ka) Pr 4 (ka) 
Pp, , ;(ka) 
+ cos kh ee 
Pr , 1 (ka) 
(5. 106) 


Finally, the total surface curren induced by the antenna and its image 


can be written as 


(0 0] (00) 
I ) = ) BP (cos@)4 >, B' P (cos ¢@) (5. 107) 
(1, Total n0 Dee =o "3 
where 
: n+] 
BL = CD 'B, (5. 108) 


Therefore, the total surface current due to both antennas is expressed 


as 


4 


[as] = 2), BP (cos @) (5. 109) 


Total n=odd 


B 


n 
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where 
= max [fn - es - coskh Pn - ie) 7 Pas 4) 2 ee RR Pus (ka) 
a) Pin. 1 tka) Pie 7 A ka) Pos y (kay Pay 1 kad 
(5. 110) 
Finally, 
I, (O)= In 3 pon : a - coskh Pn - ve : Pas 
ax n=odd Py. 4 “a P,_ 1 (ka) Pook y 
Pa + ka) 
+ coskh Dp. (Kar ] Pp, (cos 6) (5. 111) 
n+l 


5.5 Far-Field Expressions Due to a Spherical Surface Current Distribution 


The radiation vector N due to an electrical source I' is 


found to be 


N= fiftrel™ R008 v ays (5, 112) 
Vv 


where y is the angle between the radial lines to the source point and 


the observation point. 


The far-zone electromagnetic fields are obtained from 


radiation vector N: 
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4 
f 
F 
4 


arate 


aa CVT TS 


P(R, 9, ¢) 


Xx 


Fig. 5.7. Fields due to a surface current ona 
spherical ground plane. 


| 
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E, = Hh eK Ry, 


6” 42R 
(5. 34) 
_ Jou ,-jkR 
By = dar ° Ny 
Since R' =a at the surface of the sphere, the term elk Rosy in 


Eq. 5. 112 can be written as 


eikR' cosy _ ,jk a[cos @ cos 6' + sin @' sin 8 cos (¢ - 9') | 
(5, 113) 


The current density I' inside the volume integral is the 


i 
surf..ce current per unit width and flowing in 9 direction. Therefore, 


ra Van 
I' = k(0')® = k(®") [ cos @' cos OR + cos 9' sin oy - sin Q' 2] 
(5, 114) 
where x 4, and @ are unit vectors of the Cartesian coordinate system. 


The radiation vector N can be written 


7 


= 2 lan A A 
N = f k(6')[ cos @' cos ¢' x + cos @' sin g" y - sin 0'’z] 
7) 


Ot. 


eika[ cos 6' cos 0" + sin 6" sin @ cos (9 - 9')] .2 5.) g gg dg" 


oe (5, 115) 
Breaking N into three components 
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rf 


2 
N. = f J") cos @' cos 9" 
6 0 


eika[ cos @ cos 8" + sin @' sin 8 cos (9 - 9')] a? sin 0' dO’ dg’ 


Qn 
N= f 
oO 


Ose 


k(8") cos @' sin ¢' 


eikal cos @ cos 0' + sin 8" sin 0 cos(¢ - 9*)] a? sin 0' do’ dg’ 


(5. 116) 
an 
N, = -f Jk’) sino 
0 0 
elkal cos @ cos 9' + sin Q' sin 9 cos (¢ - ¢')] a? sin @' d0' as! 
Also 
2n 
J k(@')a sin @' dg’ = 1(0') (5. 76) 
re) 


Since k(@') is independent of ¢' 


n=odd 
B_ P (cos @') 
1 — __1@') ace 5. 117 
k(6") = 2na sin®' ~ “nasino’ ( ) 


Using the identity(Ref. 14 p. 407) 


gikR cosy = y (j)"(2n + 1) ) (48) P (cos w) 


n=0 


(5. 118) 
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where 


fate) =a In 4) (5. 119) 


and expanding P (cos y) into a finite series of the form 


c n 
P (cos y) = > P (cos 0) + 2 (c, COS m¢ + da sin m¢) Pi" (cos Q) 


(5. 120) 
and evaluating the coefficients c. and a by using the orthogonal 


relationship of the Legendre polynomials such as 


an 1 
m , _ 2a (n+m)! 
J J P_ (cos y) PH (cos @) cos mg sin @ d@ dg = Gnvi G@-m)t Cm 
(5. 121) 
and 
21 7 An 
J J¥,@.9) P_ (cos 0) sin@ do dg = >> [ ¥,(0, 9)] (5. 122) 
eae) g=0 
where 
¥ 0,0) = P (cos Q)cos m¢ (5, 123) 


cos y = sin @' sin ®@ cos (¢ - ¢') + cos 8 cos 9 


it can be shown that 
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an 2 
; . 4a ' ’ 
J JP, {eos vy) P, (cos 0) cos m¢ sin @ dO dg = nat Pr (cos 0’) cos m¢ 
(5, 124) 
The coefficients Cn and qa are evaluated to be 
ce = giz ml! P™(cos @') cos mo! 
m (n+m)! “n ; 
(5. 125) 
_ 9(n-m)!  m cuss ; 
i = 2 aan) P , (cos @') sin mp 


Therefore, 


P (cos v) = P (cos Q') P (cos Q) 


n 
gv (oem)! P™ (cos 9") P' (cos Q) cos m(9 - o') 


oiy (n+ m ! 
oe (5, 126) 
\' t 
ann Boos Pagay (os 9) 
= i J 2m a sin Q' cos 0" cos 9! 


ik @ cos @' + sin @ sin 0’ - 0 
eika[ cos 6 cos 6 + sin @ sin 8’ cos (0 - 0')] a? sin Q dO’ do’ 


(5, 127) 


Since 


% ‘ 
ee Wi eS a one a) (ka) P (cos 0") P (cos @) 
n=0 


n 
VY (n-m)f ym 


ee (n +m)! Pi 


(cos 0') P (cos @) cos m (0 - 0’) 
me=l 


(5, 128) 
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and 
cos m (¢ - ¢') = cos m@' cos m¢ + Sin m¢ Sin m¢' 


the contribution of the first integral with respect to ¢' 


reduces the double integral to 


’ Bo get Po 943 (CO8 9) 


" e 
i 
NG) eee gg tae) Fe) 


-1)! 4, excane 
nat P_, (cos 0") Pi (cos @) 7 - cos¢ a* cos Q' de" 


7 
= [2a B P Q' n 2n+1 
fe 2041 Po g4s (eos di nn +1) 


Fn) P. (cos 0") Pp (cos @) cos @' | cos ¢ (5, 129) 
because 


an 
f cos m¢' cos ¢' dg' 
re) 


0 for m# 1 
= 7 form = 1 


f Sin mo' cos ¢' dd' 0 for all m 


Similarly, for the evaluation of Ny , we obtain 


n ane+i 


yy 
B feos) 2 ae ay 


vv 
4 Pagel Pees 


: (5. 130) 
2 fo (ka) P’ (cos 0") P’ (cos 8) cos @' d0"| sin o 
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If we ‘et 


u nn 2n+1 
N- @ SD Bae Pag) 2 9 ate ay 


2 fn) Pr (cos Q') Pi (cos 6) cos 6" dQ" 
(5. 131) 


then 
= cos ¢ 
Ny = Nsine 


(5, 132) 


A as 
Finally, the z component of the radiation vector N can 
be written as 


a i ig k(6 ') sin 9’ eka [cos 6 cos @'+ sin 6@ sin 6' cos ¢' Ja’ sind 'do'd¢' 
0 4% 


an B P (cos 6 ') 
=--f f—Ht ee Yan +n) 4 (xa 
0 0 n= 


[P,, (cos @ ')P_(cos 6) +2 5 {n= _m): P™ (cos 6')P™ (cos 6) cos mg’ | 
fool | (n+m)! “n n 


a* sin 6'dé 'dg' 


(5. 133) 
Since 


1 
J cos m¢'d¢' = 0 
0 


for 


q 

a 
Lt 
r* 
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™ © 
= -2 ' 
2 ad d Boo 1 Pog 4 s(cv8 9 ) 


@ 
[) j*2n+1) (ka)P_(cos 6")P (cos 8) ] sin 9 "da ' 


n=0 
(5. 134) 
From the orthogonality relationship of the Legendre functions, 
1 
f P_ (cos 6) P (cos 6) (-Sin 6) d@ = 0 m # n 
9 m n 
m n 
2n+1 (5. 135) 


Tkorefore, 


22+} 


za 
! 


(oa) e 
= 2 
274 J2B 201 foo 41k2) Po 941 (C08 Q) (W3 73) (42 +3) 


ice) ° 
.\2L+1 
4a bE Bog st) Yog11 4) Pog, (C08 @) (5, 136) 
Radiation vector N ina spherical coor?inate system 
is related to the Cartesian components by 


Ng 


(N, cos ¢ + Ny sin ¢) cos 0 - N, sin 9 
(5. 137) 
N, = -N,, sing +N, cos @ 


Therefore, the transverse component of N becomes 
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N, = Neos 6 - N. sin 


(5. 138) 


Due to the symmetry of the problem with respect to z 
axis, it was anticipated that there be no ¢ component of E fields, at 
far-zone. The result shown in equatior= 5. 137 and 5. 138 confirms 
the physical phenomena. Since the direction of the E fields coincided 


with that of the N, and it was shownthat N. = 0, only the @ com- 


g 
ponent of E exists. 


Finally, 
: -jkR 
_ jk, ] 
a ~  ~GaR t 
(5. 139) 
; -jkR 
jkn,e J Ze 
aR [Ncos Q- N, sin 0] 0 


Equation 5. 139 is the expression for the far-zone electric field E. 


due to the spnerical surface current that was; in turn, induced 


by the monopole and its image. 
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5.5.1 Evaluation of 1. From Eq. 5. 131, R is defined as 


1" 
: wn 2n+1 ‘ 

. J PD Boost Pog 1(608 9) ey } n(n+l) 4, ka) 

Pi(cos 6") P (cos 6) cos 9' dé' (5. 131) 


since 


qa™ P (x) 


ax™ 


£m) = (-)™ (1-2)? 


1 d P (cos 6") 
P (cos 6 ) de" = “acos a") d(cos 6 ) ’ 


if we let x = cosv', the integral can be rewritten as 


TN = 2 J Zs Boo a1 Pops a i" aay 4) 


d P (x) 


x dx (5. 140) 


1 
P cos 6) 


The evaluation of this integral depends upon the evaluation 


of the integral I 


1 d P_(x) 
I = I x Pop , s(x) 


(5. 141) 


From Ref. 17 
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d P (x) 
oo = (2n- 1) P 1) + (2n-5) P3(*) + (2n-9) Pid 5(x) ton, 
(5. 142) 
also 
(2n- 1) x Pe 1) =n P(x) + (n-1) P. 9(X) (5. 143) 


From (5. 142) and (5. 143) 


dP _(x) 
x + = [n P (x) + (n- 1) P_o(X)] + [(n-2) P -9(X) + (n-3) as 4g] 
+ [(n-4) Pes 4) + (n-5) P 5(x)! 4 ie ou 
= MP (x) + (2n-3) P_o(x) + (2n-7) Pg) + (2n-11) Pox) +. 
<n/2 
= nP (x) 4 L (2n-4q+1) P aq(*) (5. 144) 
Thus 
1 
is J Pel) [n P (x) +(2n-3) P(x) + (2n-7) Pcs) ties Jo 
(5. 145) 
From the orthogonal relationship of the legendre polynomials, 
J 
3 P(x) P (x) dx = 0 for mn 
(5. 135) 


= 2 for me=n 
~ Bn+i etc: 


159 


Using the orthogonality shown above, the integral I can be 


written as 
i 1 0/2 
I= io nPo, 10) Pytad a+ f 2 (2n- 4q+1) Pop 10) P.,_ 9 (x) ax 
- - q=1 
26 
n-(20+1 n- 2-1 2 
= (2€+1) et + (2n- 4 pea Tecan 1) Bes + 1 
(20425 8 Sn/2 
n-(22+1) 
2 $53 +2 2 6 2q- (2041) (5. 146) 
where 
ie = 1 for m=0 
0 for m#0 
Then, 


os 2041 2(2041) +1 YQ 1 2( 2 +1) 
Y= 2a d Boosi) GESTED) Foest*) Pop 10089) “grag 


n 2n+1 
+ 4a py Boo 41 J nine) antl 4 (ka) P * (cos 6) 


li 
bo 
fo 
ub-48 
iss] 


1 1 
9¢ 416) Tri 4 op 4108) Po, (C08 6) 


’ on 2n+l 
+ 4a 1H Boo j PRE crs) Rat) 4 a ) pi 7 (C08 9) ; (5. 147) 
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Also from the orthogonality condition of the second integral in I; i.e., 


n-2%-1 
a 
it follows that 
n= 2£+2q+1 
where 
£ = 0, 1, 2, 
q = 1, 2, 
n = 2843, 2£+5, me: 
Finally, 


e bag (ka) 
1 - Qa y B 2041 qoos1 pl 


- 2041 S57 Pay 41 (e084) 
£=0 
© (ve) 2 
.2m+i _ (4m+3) 1 
+ Po Bon 41 se J (2m+1\(m+1) tomsi®) Po my 108 9) 


(5. 148) 


The far-zone electric field E 9? due to the surface current on the con- 


ducting sphere induced by two monopoles, is: 
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‘ene ikR 
~ik Ip © 4 28+1 fon 41) 1 
Xo, aR [acoso( 5 Pats) “EST Paps s(0084) 


+ 


Y « .2m+1 (4m+3) { 1 
ds Patel ae (2m+1)(m+1) fomistk®) Pamss(o5) 


+ 


F 2041 4 
2a sing oy Boost j Gop 1088) Pop , (608 | 


which becomes 


-jkR 
Ny © 


: B, ; 
9° aE ae lb, 4,(sa) Pi(cos 6) - = 4 (ka) P4(cos 6) 


: “ 4 ska P5(cos 6) | 

+ Bags! , (e 4 ka) P3(cos) +45 4.) P5(cos 6) 
- 5 4 (xa) Pr(cos 8) + «0... \ 
+ B, (75 4 a P5(cos 6) - x 4) P7(cos 6) + mo 


+ Be (- z Fox) Pr(cos 6) +... ) + kasin@ rac P (cos 6) 


‘ ‘ \ 
és Bs Ta(ka) P,(cos 6) + B, t stra) P,(cos Q)-.. | | (5. 149) 


The ratio of the first term to second term is, in general, 


much greater than 1 as was shown in Tables 5. 1 and 5. 2 so the series 


converges rapidly. To simplify the expression, an initial term is 
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taken to approximate the infinite series. It will be shown later than the 
subsequent terms do not contribute significantly toward the exact value 


of the series. 


The simplizied form of E Q is 
S 
Loe Ue) -B resting Bi cos 6 sind ; ; 
(5. 150) 
where 
2 B. 
BLY 
max 


The numerical calculation of the far-zone fields E 6 are 
Pe) 


performed with the above equation (5. 150). 

5.5.2 Numerical Evaluation of a Radiation Pattern. Since, 
in Section 5.3, it was shown that the far-zone electric field due to two 
monopoles separated by d = 2a + h, measured between center to center, 
has an expression as shown in equation (5. 69); i.e., 

; inlme on 


“64 7 oak 


cos(kd cos 6) - cos k(d-a) cos(ka cos 9) + cos 6 sink(d-a) sin(ka cos 6) 
Sin 6 


(5. 69) 


The total electric field E, at far-zone is 


E, = E, +E, (5. 151) 
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This theoretical evaluation of a far-zone electric field due 
to the monopole above a semi-spherical ground plane on an infinite 
ground plane was conducted mainly to evaluate the effect of the detailed 
size of the semi-spherical ground plane, although its radius was much 
smaller thana wavelength. Here, a majority of experimental studies 
were performed with models that have the dimensions of a, corres- 
ponding to the radius of the sphere, less than or equal to a quarter of 
a wavelength of interest. 

The expression E 9 is in the form of an infinite series 
where the coefficient BL and n(K2) are constants for a given geo- 
metrical size of the ground plane and a given frequency. BH is a com- 
plex function of the radius of the sphere and the antenna height in terms 


of wavelength; i.e., ka and kh. 


When the argument of a spherical Bessel function is much 


smaller than 1, it can be evaluated using the formula 


‘ n 
~~ Z ~ 
Fn(2) © aye for = 0,1, 2... 
z<<1 (5. 152) 


Although the above formula is for z<<1, the ratio of these functions 
hetween n=1, and n=3 or n=5 are fairly close to the actual ratio 
of these functions evaluated with exact values of fa?) , 


Since the far-zone electric field EQ due to the spherical 


surface current is in the form of an infinite series with only odd n's , 
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the ratio of the successive terms and the cature of its convergence is impor- 
tant to justify taking only the first few terms in the actual evaluation 
of the field and the radiation resistance. 
Ratios of the f (x) to 400 and 4,(m to 4) are 
evaluated by using the approximate formula (5. 152), and by using the 
exact values for x = 0.5, }. and i.5. Results are tabulated in Table 5. 1. 


Because the ground plane sizes used in the study were 


a 
2 


spherical ground plane equal to or less than one half of a wavelength, 


limitedto ka < which corresponds to the diameter of a semi- 


the xatios only of those spherical Bessel functions with an argument less 


than or equal to 3 are of interest to us. 


able 5.1. Ratios of 4. ' 
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Further study is necessar; to determine the magnitude of 
each term of the series. Each term, in addition to 4,00 , contains 
Bo the constants evaluated for each given size of the ground plane 
radius and the antenna height. Ratios of By to By and By to Be can 


be obtained from Table 2 for several typical values of ka and kh. 


[a no [aes at [ant won 


1, 298 - j0. 742 0.808 - {0.078 0.342-j as 
B 247 +50. 738 -0.057 +j0.078 0.052+j 1.30x 10° 


-0.046+) 4. 28x 107° 0. 0287+ j9. 66x 107° 


0.031+j1.11x 10° 


Table 5.2. Coefficients B," 


Finally, the be havior of P cos Q) and P ‘(cos 6) are 
graphically shown in Figs. 5.8 and5.9 for O<@< 90°. From 
Fig. 5.8, the magnitude of P(x) is to be always less than or equal 
to 1, and the ratios of P *(cos ®) for n=1 and 3 or n=3 and 5 
are always less than 6. 

hk can be concluded, with the considerations given above, 
that the ratio of the successive terms of the series between n= 1 


and n=3 is always lessthan 0.1. Therefore, it is jusitifed using 


only the initial term from each infinite: series in Eq. (5. 149) or (5. 153) 
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Fig. 5.8. P (cos 0) versus Q. 


P.. (cos @) 
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Q@( Degrees) 


Fig. 5,9. Plicos 9) versus @. 
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co 4 (ka) 
. 2441 2k+1 

p J Boral £+1 ee (cos @) 
i 4= 

vB 2 j2m+1 (4m +3) bixa) po 

Paes! (2m +1)(m+1 foms 2m+1 (cos 9) 
£=0 =Q+1 
(5. 153) 

= -2641 7 

p Bo git : Foes) Po 941 (€°8 9) 
: Thus, 
L ' -jkR 
inl? 
4 ] anR 
i 


sees - cosk(d-a) cos(kacosé) + cos6sink(d-a) sufhacoe 
Sin 6 


o 


((ka) B,"cosé sing 
: i| Nees: iE 4 : 4 (5cos® 9-1) sta) 
(5, 154) 


The plot of the magnitude of the term inside the bracket is 
a function of 6 and is a radiation pattern for a particular a and h. 

A few of the patterns are calculated and are plotted below in 
Fig. 5.10 in order to compare with the experimental results. 

Figure 5. 10 displays a set of theoretical radiation patterns 
E 9 for a monopole and for its image in the upper hemisphere. This is 
the plot of Eq. 5.154. Since the radiation patterns are symmetrical 


along the 9 = 0° axis, only the right half of the complete radiation 


boa aca. | 


FORT TEP TOT THR 


d={.25m 
EX 
Q=2.50m . . 
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Q=3,125m fs ; 
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eee ee 
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@-0" 
'p 
h-2.5m : ; 
h=A/4 ot 30MH2 a 


Fig. 5.10. Theoretical radiation patterns |E 9 |? 


; 
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pattern on the x - z plane is plotted in the polar coordinates. 

A minor lobe off 9 = 0° axis increases with the size of 
ground diameter. The changes in the ratio of the intensity of the major 
to minor lobe is shown in this figure. 

There are a total of 9 patterns plotted in Fig. 5.10. Each 
row is a plut with a constant distance of separation, and each column 
with a constant frequency. The monopole antenna was fixed at a quarter 
wavelength for the highest frequency in the band. Each pattern contains 
infor mation on the antenna length and the separation distance between 
the two antennas. 

Comparison between the patterns theoretically obtained 
for a monopole and its image separated by a diameter of the spherical 
ground plane, but without the presence of the ground plane, to those 
with the ground plane shows that adding the ground plane moved the 
null position between the major and minor lobes toward 6 = 90° axis 
at higher frequencies. Also at lower frequencies, the beam width was 
made to become broader by adding the ground pvlane. Comparing the 
theoretical results to the experime al results also show that the 
general behaviours of radiation patterns are similar even though actual 
amplitude do not seem to match exactly. ‘his can be explained by 
pointing out the fact that the initial assumptions made concerning the 
current distributions on the antenna do not auite agree with the actual 
current distributions. Also a finite ground plane used to simulate an 
infinite ground in the scale models caused the whole pattern to be shifted 


toward 9 = 0° exis. These facts are theoretically shown by Leitner and 


Spence (Ref. 3) in their study on a monopole with a finite ground plane. 


ae ae Bn SS Oe ea NE ERE a a eR ON ME ie Win, AD Din ret Le 
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5.6 Radiation Resistance 

In this section, the radiation resistance is calculated as 
a function of the antenna length and grourd plane diameter. The radi- 
ated power is calculated using the far-zone field expression that was 
obtained through an approximation of the infinite series. The far- 
zone field Ey was given in Eq.5, 154. 

The convergence of Eq.5.153 for ka<7, justifies 


using only the first few terms to obtain E 0° The far-zone electric 
Ss 


field due to a monopole of length h and a hemispherical ground plane 


of radius a is 


ine KR 
Bi, oes a 
is) 27R 
i cos(kd cos @) - cos k(d ~ a) cos(ka cos 0) 
m sin 0 


+ cos § sin k(d - a) sin(ka cos Q) 
Sin @ 


(ka) Bi cos 9 sin @ 
- j 5 


f,«) -46 cos? @- 1) 4,0) (5. 184) 


The first and second terms of the equation in the brackets 
are due to contribitions from the currents on the two monopoles and a 


spherical ground plane, respectively. 


a ee Se NT eg ie, ae Te Rae ee a ae ee ae aie ae Ne SSS eR Eee ee 
6 SO er eR Nope tT OM Uy we Peg 2 Re, a SN Eg ee ae ee ee ee Ne To a A 
ee RR Eg Ee ee ee ee a a ey 
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F.(0) = cos(kd cos @) - cos k(d - a) cos(ka cos @) 
ES A sino 


i Se ae ee aaa ke hae ee 2h tit d ts Ba) 
- 


+ cos 9 sin k(d - a) sin(ka cos Q) 


sin 0 (5. 71) 
and 
o q 3 e 
(ka) cos @ sin @ | d (ka) +z (5 cos” @- 1) data (5. 155) 
F4(6) PE oe ae a re a ee ee 
2 
Then, the Eq. 5. 154 can be written as 
“ -jkR 
My & ; or 


Letting By =a,+ jp 1? Eg can be written as 


ine? 
0 ; 
Bo = aR I [Fi ar Oar F,0)| 


jn e JKR 
en, oa a ee 5. 157 
oa In | [Fy By #516) Jay F,(0)| Sey 


The far-zone magnetic field intensity is, from H = ia (Rx E) , 
fo) 


. = jke 
je 
H. = 


¢ 97R ey [Fo + BFA) | - j ay F,(0)| (5. 158) 


Ne Kiara IR a OE, a EE RE TN ES SE RN Fee ORE OPS Te SETS re se agen a a Aa a aka «| 
ae ee St te ° 
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. * = j elkR : 
; Hy Sie Ts [Fo By F4(0)]+ ja, F,(0)| (5. 159) 


Radiation resistance R. is defined as 


2 


n I 
1 om 2 2 2 2: 
7 g dn? R? [[F,@ +By F4(0)| +a," Fy (@]R sin @ dO dg 


7 
= ras § F 1° (0) + 28, F,() F@) + 8,7 F*@) +a,? F,* 


Sin @ dO dg 


1 
= 204 2 2 2 2 2 | < 
7 SF, (@)+ 26, F,(@) FL(@) + (2,7 + 2,7) Fy? (@)| sin 0 ao 

(5. 160) 
Since No = 120 7 in free-space, 


T T 
7 2 : : 
R = 60 J [F, () | sin @ dQ + 120.8 f F ,(@) F,(8) sin 0 ao 


T 
2 2 cr 24 : 
+ 60(a, + By J Fy {Q) sin © d@ 
(5. 161) 


R,_ =R, +R, +R (5. 162) 
r k 


where 
‘ 2 
R, = 60 f F,7(@) sin @ ao 
1 ) 
tT 
= 5. 163 
R,, 120 8, J F ,(0) F,(@) sin @ de (5. 163) 
2 y 2 
_ 2 ; 
R,, = 60(a,7 +8, ) SF (8) sin @ do 


5.6.1 Evaluation of R.- From Eq. F (0) is given as 
1 


cos(kd cos @) - cos k(d - a) :os(ka cos @) 
EF) = ee eee 
1 sin 0 


+ cos 9 sin k(d - a) sin(ka cos Q) 
sin 9 
(5. 71) 


Therefore, 


7 Aas A PG LTE RTT RT EEO TF TS i. 
me gee pee eG ERORTERLAT OEE EO INTE ET PIT 2 er alias 

OPT SSE BETO SE SEC PD 

FLT FTE 

poe 
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Ee ‘Grae (kd cos 8) + cos? k(d - a) cos* (ka cos @) 
sin 


+ cos? @ sin? k(d - a) sin? (ka cos 9) 
- 2 cos k(d ~ a) cos(kd cos 9) cos(ka cos @) 


+ 2 cos 9 sin k(d - a) sin(ka cos Q@) 


2 cos 9 sin k(d - a) cos k(d - a) cos(ka cos @) sin(ka cos 0] (5. 164) 


1 
is 2 : 
R as 60 J [F, @)| sin @ d@ 


<3 ae 


sino 


6 e f gos lk ? (kd cos 9) cos 9) a0 
6 


(2) = cos* Kd a) f S08 boss We 208-9) cg 


T 2 ee 
a. Sesela cos" @ sin“ (ka cos @) dd 
= sin” k(d - a) J ~——snd 


1 nr 
() & soreowkid a) f cos(kd cos S} coaika cos Q) d@ 
6) 


() = 2-mx(d-a) je Soa dco 0) silat tvs 0) ae 


O 


ee e 
@ = -2 cos k(d - a) sin k(d - a) fees cos @ cos Sika cos 8) sin(ka cos @) dQ 
Oo 


176 


Let uw = cosé, dp = -sin@ dé and d = a+h and 
solve integrals from @) through © separately. Then sum of the results 
will give ae . In evaluating each integral, logarithmic singularities 
appeay These singularities will later cancel. However, in order to 
show the convergence of the integrals, the limits of the integration have 


been changed from yp = -1 and pw = 1 to wp = -l+e and w= l-e 


and then € was made vanishingly sma!l. 


Then, 
Ive cos”? k(a + h) uw 
(1) = lim { aT ean Gelb 
€~-O0 ~}+e 1- pu 
Since 
1 _1 1 1 
- Mh 2\1l+p 1-4 

and 


\ 
cos*k(a +h) p = 3 (1+ cos 2k(a + h) H) 


€~0 - 1l+e 


ete ie Oe ee oe 1 +cos 2k(a +h) p)d 
G@)= qlim J i+p l-p sd 


b b 
F F 


Therefore, 


1 
hostels 1 + cos 2k(a + h)py 
(Qa) = 9 lim if a. Pg du 


1 1 : 
= ; hn fi rear + f achliadt “Sabin h) 4 au| 
E~OL- l+e - l+e 
After the change of variables, 1+ 4 =u for the first integral 
and 1+u=u and 2k(a+h)u =v for the second integral, @) can be 


written as 


2 4k(ath; 
@) = + lim f au | cos 2k(a + h) f saa dv 
u v 
e~7OLle 2k(a+h)e 


(n2- lim (ne + cos 2k(a + h) 


4k(ath) |. 7 
+ sin 2k(a + h) i ~~ | = ; 
me) 


2k(ath)e 


le fix + h))- lim Ci (ox1a + ne) + Sin 2k(a + h) Si (x(a + h)) | 


€~-0O | 


where 


X,- 2 
Si(x) = f a dv = sine integral 
Oo 


x 
Ci(x) = f ——s V dv = cosine integral 

eo) 
In the solution for (1) , Singularities associated with logarithm 


and cosine integral are left as they are and it will be shown later that they 


cancel. 


Using a similar technique, it can be shown that 


l-€ a 
(2) = cos* (kh) lim i: cae es du 
emo -ire =o - 


d 
2 


cos? (kh) ) én 2- lim {ne + cos 2ka [oi ka) - lim Ci(2ka «| 
6-0 €—0 


+ sin 2 ka Si(4 wat 


RNS FEES I BETES TE BER TS FPS 


© 
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i 
(sin? kh) tim) f—-#2_ Sin? (ka w) ds 
€—0 e 1 - ue? 


© 


[+e 


adding and subtracting sin? (ka 4) in the numerator. 


l-e 


. 2 2 oj 
= (sin? kh) lim f Sit (kau) (= ") sin? kaw gy, 
EO = Ire l- yu 


l-€ - 2 I-e 
= (sin? kh) lim f Sin ka yt du -f sin* ka pb au| 
e-Ol-lte l-u - l+e 


The result of the integration is 


@) = 5 (sin? wy | in 2- lim ene -2 (122i) 
€—0 eo 


\ 


- cos 2ka [ew ka) - lim Ci(2 ka | - Sin 2ka Si‘4 ka) 
€~0 


l~e 
~2cos(kh) lim f $08 K@ +h) u cos(ka y) 


2 du 
€-+0 -]+e¢ 1- i 


uSing an identity that 


cos AcosB = a [costa + B) + cos (A - B)| 


- cos(kh) lim 
€—0 


i cos k(2a+h)u a, i cos {kh 
l+4u u 1l+yu 


- ite 


{ 
- cos(kh) } cos k(2a + h) Bec + n)- lim Ci (sea +h) I 
6 


+ sin k(2a +h) Si (2K(2a + h) + cos kh eve kh) - lim Ci(kh € 
€—O 


+ sin kh Si(2 kh) 


1- ’ 
2 sin kh im i u. cos k(a +h) uw sin ka p du 
€~0 -l+e i- p? 


using the identities 


and cos A sinB = [sina + B) - sin(A - B)| 


1 
2 


and also po du = Le F(u) du when F(u) = -F(-p) 
, i- a ite 


the result of the integral G) becomes 


= sin kh} sin k(2a +h) oi( ax + n) - lim Ci (k(2a +h) J} 
€+0 
7 
- cos k(2a + h) Si (2x20 + 1) - sin kh cca - lim Ci(kh a 
E-0 


+ cos kh Si(2kh) 
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Similarly, 


l-e ; 
© = -2coskhsinkh lim f HcoSkausinkay 4, 
EO -1+e L- p? 


2 + sin 2 kh { cos 2 ka Si(4 ka) - sin 2 ka lows ka) - lim Ci(kh ol 
€—0 


Cosine integral Ci(x) can be expanded into a power series 


suchas (Ref. 18, p. 232). 


2 42 gen ; 
Ci(z) = y+ fn Z+ > “Sn(an)t (5. 165) 


If argument of the cosine integral, z, is small 


lim Ci(z)= y+£nz (5. 166) 
Z—0 


where 


y = 0.577 = Euler's constant 


Therefore, lim Ci(2k ae), for instance, can be broken up into two terms; 
rome) 
ye =a 


lim Ci(2 kae)=y+ €n(2ka)+lim  Ln(e) (5, 167) 
670 €—0 


a 
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Applying the technique shown in Eq.5.166and adding the results @) 


through (6), it can be shown (Appendix E) that 


R 2y + 2en(2 kh) - 2 Ci(2 kh) - 2 sin® kh (1 - 


a 


30 


sin 2 ka ) 
2 ka 


+ cos 2k(a + h) [ve ka) + ci(' k(a + 1) +2 £n k(2a +h) - £n(2 ka) 


2si(2¥4 +») 


If we imzase a condition that a=0 and h=m ; 


gn 2k(a +h) - 2Ci [ax(2a r »]| + sin 2k(a+ h) su ka) + si (tk (a+ ) 


(5. 168) 


where 


m is an odd integer, the result matches with that of Stratton's (Ref. 14, 


p. 444) where he found a radiation resistance of a linear dipole at reson- 


ance. 


When a=0O andh= + , Eq. 5. 168 reduces to 


R.. = 30 [y + £n 27 - Ci(2n)| 
1 


which is the exact replica of Stratton's expression 
R = 30[gn2ny*- Ci ama] with m= 1 


y*= 1.7811 


It can further be shown that when h = 0, R.. = 
1 


(5. 169) 


(5. 170) 


0. This, 


of course, must be true for physical reasons and this gives andther 


confirmation of the correctness of this theoretical result. 
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5.6.2 Evaluation of R 
T? 


From Eq. 5. 163, 


ee) 
ul 


T 
1208, { F,(@) F,(0) sin 0 do 
oO 


E cos(kd cos @) - cos k(d - a) cos(ka cos @) 
120 g f ee Ea 
1 sin ® 


Oo 


, £08 8 sin k(d - a) sin(ka cos 9) | 
sin 0 J 


(5. 171) 


Let 


(5. 172) 


Then 


F,(9) = G(ka) sin 9 cos @ + H(ka) sin 9 cos® @ 


R can ow be expressed as 
2 


r 


BITE gE OEP Ee ET ne 
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R. = 120 A, f cos kd x - cos k(d - a) cos ka x + x Sin k(d - a) sin(ka 2] 
2 


[ Gea) x + H(ka) x*] dx 


Every term in the above integral is an odd function of x with 


respect tox=0. Therefore, the result of integration becomes zero. 


therefore 


5. 174 
Ty ( ) 


5.6.3 Evaluationof R. 
T3 


a aa 


T 
= 2 2 2 5 
R. = 60(a, + B,7) J F,~ (8) sin @ do 


1 2 
60(a,* + B,*) i [ Gta) sin @ cos 0 + H(ka) sin @ cos*o | sin @ dQ 
o 


1 
= 60(a,* rm B,*) f [ G (ka) sin? @ cos?@ + 2G(ka) H(ka) sin?@ cos‘@ 
0 


+ H* (ka) sin? @ cos ‘0 | sin 0 dQ (5, 178) 


Letting 


cos@ = x -sin@ dQ = dx 


sin*@ = 1 - x? 


1 
60(a,7 + 8,7) f [G* (ka)(d - x*) x* + 2G(ka) H(ka)(1 - x?) x* 
-1 


+ H® (ka)(1 - x*) x*] dx 


1 
60(a,* + B,*) J G? (ka) x? + [ 2c(ka) H(ka) - G? (ka) | x* 


+ [ ne (ka) - 2G(ka) H(ka) x® - H® (ka) x® 


60(a,7 + B,*) G? (ka) (F) +[2c(a) H(ka) - G? (ka) (2) 


+ [2 (ka) - 2G(ka) H(ka)|2; - H? (ka) ) 


60(a +8 a Te G? (ka) + a G(ka) H(ka) + = H? (ka) 


240(a,° +B re G? oe. » Sate) Hike) r ml) 


(5. 176) 
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Finally, the total radiation resistance of the system with a 
monopole and a hemispherical ground plane is ubtained by summing 
,»R and R_ . However, this sum represents a radiation 
Ty To rT, 
resistance of the actual antenna system and its image. Because only 
a half of the power evaluated previously is actually radiated, mono- 


pole and hemispherical ground plane radiating into a half Space, the 


final radiation resistance is 


4 
Rad. Resistance = ZR, + Rr + R.) 


15 2, + 2 £n(2 kh) - 2 Ci(2kh) ~ 2 sin? kh (2 - Te ma) 


+ COS 2k(a + h) [eve ka) + ci( Ata + 1) +2 ¢n k(2a +h) 


~ £n(2 ka) ~ en Qk(a + h) - 2 Ci [ax(2a + »)] 


+ sin 2k(a +h) | ka) + Si (+x(a + n) - 2 Si (24(a + ll 
I 


2 2,|G*(ka) | 2G(ka) H(ka) (ka) 
ds a | Ae te ~ 3 


(5. 177) 


From Eq. 5.177, radiation resistances of several different 


combinations of ka and kh are numerically caleuiated and tabulated. 


The result of numerical evaluation of Eq. 5.177 is given in 


a ,raphical form in Fig. 5.11, where radiation resistance is plotted as 
a function of frequency. Each curve represents a different size of semi- 
spherical ground plane. When the radius of a ground plane is zero, 

the resulting radiation resistance corresponds to that of a monopole on 
an infinite ground plane. 

Results shown in Fig. 5.11 indicates that there exists 
definitely a peaking effect on the radiation resistance as the radius 
of the ground plane is changed. Comparing with the monopole resonating 
at 30 MHz over an infinite ground plane the radiation resistance becomes 
larger with a semi-spherical ground plane well below the resonant 
frequency. The peaking seems to occur at ka = 1 and this conclusion 
has been drawn mainly from the results of numericai calculations. 

Fig. 5.12 and 5. 13 compare these theoretical results with 
experimentally obtained input resistances. Of course, we are not 
comparing the same resistances, namely the radiation resistances. 
However, assuming that the loss is small, the input resistance should 
be similar to the radiation resistances. Some of the discrepancies 
shown in this comparison can also be explained with the discrepancies 
in the assumed current and the actual current on the antenna. However, 
the peaking effect is shown to exist using a smail semi-spherical ground 


plane with a monopole. 
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Fig. 5.11. Theoretical radiation resistances for various 
values of ground plane size 
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Fig. 5.12. Theoretical radiation resistance and experimental 
input resistance for a monopole with a hemispherical 
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Fig. 5,18, Theoretical radiation resistance and experimental 
input resistance for a monopole with a hemispherical 
ground plane 
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Since radiation resistance is defined as the ratio of radiated 
power to a square of the input current at the terminal where it is mea- 
sured, the current must be that current at the terminal and resulting 
,adiation resistance will be the radiation resistance at that point. 

At the beginning of the theoretical analysis for this problem, 
it was assumed that the current distribution on the antenna itself is in 
the sinusoidal form. 


From Eq. 5. 59 , the current on the antenna is given as 


I(z) = Dt sin k(d - z) 


where at Z:a, wink(d- z)=1 at 30 MHz. 


Therefure, the input currents at the base of the antenna 
(z = a) would be differen’ from I ay for frequencies lower than 30 MHz. 
We must use these curvents to evaluate the radiation resistance at the 
input terminal for each :reqvency. 

It seems as thuugh one can have a radiation resistance as 
large as desired by making the feedpoint located at a point where current 
is approximately zero. One examule weuld be a full wave dipole which 
has a zero input current at the center of the dipole. Theoretically, 
then, it has an infinitely large rac iation resistance. However, in 
practice an actual antenna is net iafinitesimally thin, which was the 


assumption used for deriving a «eoretical result, and the current at 


a minimum point is not zero. Nevertheless, the radiation resistance 
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at a current minimum may in practice be very large, possibly 
thousands of ohms. 

In view of the above argument, the radiation resistance 
cannot be used as a sole measure of an antenna performace. The 
reactive component of an input impedance is another factor which should 
be considered. When an input current is fixed at a certain value, both 
the resistive and reactive components of the input impedance determine 
the voltage required to maintain the current. Therefore, it is necessary 
to consider a ratio of reactive component and resistive components, 
or the Q of the antenna in order to decide whether increases in radiation 
resistance actually increase the antenna efficiency. Input current and 
voltage determines the input power and radiation resistance determines 
the radiated power. The ratio of these quantities determines the efficiency 
of the antenna. 

In addition to the efficiency considerations, the Q of the 
antenna system determines broadband tuning possibilities of the antenna. 
When the antenna is electrically short, in particular, the reactive 
component is usually very large and the resistive component is very 
small giving very high Q values. When the reactive component is large, 
a wideband tuning becomes very difficult, because one has to finda 
way for a wide range variation in the conjugate reactive component to 
cancel out the large radiative reactances. Also, a small value of resistive 


component will cause a difficulty in transforming up to a characteristic 


impedance of transmission iines normally available. 


The ratio of a capacitive reactance to a resistance of the 


antenna is defined as Q, for a monopole above an infinite ground 
plane, Q of the antenna increases rapidly as frequencies go down 
below resonance. Typically, for a monopole one quarter wavelength 
long at 30 MHz, with the antenna or diameter ratio of approximately 
400, Q at 5MHz becomes 1400 and gradually reduces to 0 at 30MHz. 
The resistive component increases from 2 ohms at 5 MHz to 36 ohms 
at 30Mhz. However, from the experimentally obtained data, Q of 

a monopole one quarter waveiength long at 30 MHz placed above 

a finite circular disc ground plane of 2.5 meters in diameter which 

is placed 1.25 meters above an infinite natural ground is shown to 
vary between 43 to 0. It is mainly due to smaller reactive components 
and large resistive components over the frequency range of SMHz 
through 30MHz for a monopole over a finite ground plane. 

The resistive component of the input impedance contains 
loss due to several sources in addition to radiation resistance. In 
order to make an educated judgment on the variation of the resistive 
component as a function of frequency whether it can be attributed to 
an increase in radiation resistance or not, a theoretical result of 
a radiation resistance for a mon ole over semr spherical ground 
planes has been obtained. Using these theoretical results of radiation, 
resistance and the experimentally obtained reactive components, Q 


of the antenna has been obtained between 5MHz and 30MHz. It shows 
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that the Q of the antenna varies between 25 to 1 with Xx. at 
10MHz of 850 ohms compared with 900 ohms for a monopole above 
an infinite ground plane. 
A conclusion can be drawn from the above studies that 
by using a monopole with a finite ground plane above an infinitely 
large ground plane, a Significant advantage can be realized in rad- 
iation resistance as well as a broad band matching at the input terminals. 
Through this study, a qualitative conclusion of an antenna behavior 
with two ground planes has been possible. However, an exact 
behavior of a monopole with a disc ground plane above an infinite 


ground can only be predicted with a new mathematical model. 


fT 


CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 


MSIF RAT aE TR PFI 


6. 1 Conclusion 
One of the major objectives of this study was to determine 


how a finite ground plane used in addition to an infinite lossly ground 
affects performance of a monopole antenna. Attention was mainly 
directed toward the effect of ground planes equal to, or smaller than, 
a half wavelength in diameter. If there is tobe any advantage in 


uSing these finite ground planes, they cannot be physically bulky. At 


| 


frequencies 5MHz to 30 MHz, where the study was performed, a half 
wavelength at the upper end of the frequency band is equal to 5 meters. 
In studying this antenna system, the length of monopole 
was fi.ed at 2.5 meters, which is a quarter wavelength at 30 MHz. 
Impedance measurements show that the real part is a strong function 
of both the ground plane diameter and its location. Also it shows a 
peaking effect at particular frequencies, which depends upon the 
diameter of the ground plane and its location above the infinite ground. 
Using a theoretical model with hemispherical ground 
planes sunstituted for flat discs, it was established that the sharp 
increase in resistive component was largely due to an increase in radia- 


tion resistance. From the theory, it seems that at least one peak 


exists at frequencies where ka -1. It was not possible to draw a 
clear conclusion of this nature for an antenna with a finite disc ground 
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plane. When ka = 1, the distance from the base of the antenna to an 
infinite ground a ong the path of the spherical surface is a quarter 
wavelength. This provides a voltage maximum or current minimum 
point at the base of the antenna, regarding the infinite ground plane 

as a zero potential surface. Therefore, a larger input impedance or 
larger resistive component is realized at ka =1. However, based 

on the experimental results, a general conclusion that a peaking effect 
occurs at frequencies approximately determined by (P + a): C where 
C is a function of the ratio ea. It can also be concluded that the Q 
of the antenna is generally muh lower than that of the same antenna 
over an infinite ground plane, providing some advantage in designing 
a tuning network at frequencies well below resonance. 

Curves showing the input reactance with a finite disc- 
ground plane of various sizes exhibit a region where the reactance 
variation as a function of frequency is remarkably small. This 
phenomenon is usually associated with a resistive maximum. Clearly, 
this type of frequency response is very much better, from the point 
of view of broad-band operation, than the response of either « half- 
wave or a full-wave dipole. With an appropriately designed reactive 
matching network the standing wave ratio on a line terminated in 
the antenna with a finite ground plane over an infinite ground can be 
minimized over a wide frequency range. 


Because of the assumption made in the theoretical studies 


that the current distribution on the monopole is in sinusoidal form 
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and independent of the surface current on the ground plane, actual 
current measurements were performed using scale models. Current 
distribution at the higher end of the frequency band approximated 
closely the initial assumption of a sinusoidal distribution. However, 
as frequency is decreased, particularly for the monopole with a hemi- 
spherical ground, the actual current deviates appreciably from the 
original assumption. At the lower end of the frequency band, the 
current measurement shows almost a constant amplitude similar to 
Hertzian dipole. 

The theoretical results obtained here for the radiation 
resistance assumed that the current distribution had a triangular 
form at low frequencies. These results would be improved by 
assumption, in the theory, of a current distrivution more closely 
approximating that observed experimentally. This effect is much 
more noticeable at the lower end of frequency band. Experimentally 
measured input resistances confirm these facts. In Chapter 5, a 
comparison is made between the theoretical and experimental radia- 
tion and input resistances. The experimental values are always 
larger than the theoretical values. 

From the results of a theoretical study using a hemi- 
spherical ground plane and the experimental results obtained for a 
monopole with a finite disc ground plane above an infinite natural 


ground, it can be concluded that a marxed increase in radiation 


resistance results with a finite disc ground plane, below the resonant 
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frequency. 

The major effect of the finite ground plane upon the antenna 
radiation pattern is an emergence of side lobes when the ground plane 
diameter becomes appreciable compared with a wavelength. In general, 
the appearance of a side lobe is not desirable, because the 
power radiated is actually deviated from the main lobe where it should 
be concentrated. Also the beam width measured between 3db points 
from the position of the maximum amplitude is generally broader 
than that of an antenna on an infinite ground plane. Theoretically, it 
has been found that a spherical surface current makes only a very 
Slight contribution to the far-zone electromagnetic field. 

In summarizing the results of this study, it is concluded 
that definitely improved performance can be obtained with an electri- 
cally-short monopole antenna by operating it with a physically small 
finite ground plane above natural ground. This improvement occurs 
in the radiation resistance and with respect to matching network 
design for a broad-band operation. These results have bearing on 
the design of short monopole antennas for operation of ground- 
based vehicles, in which case the antenna may be considered to be 


mounted on a small ground plane above natural ground. 


6.2 Recommendation for Future Work 


First, in the analysis of this problem, the infinite ground 


below the finite hemispherical ground was assumed to be perfectly 
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conducting, permitting use of a conventional image technique for 
solving the problem. The loss due to the finite conductivity of the 
natural ground has also been neglected. However, in order to isolate 
an exact radiation resistance from the measured impedance shown 

in Chapter 2, it is necessary to take the loss due to the natural ground 
into account. 

Also, a new theoretical model consisting of an antenna on 
a small ground plane above a larger ground plane, could possibly give 
more accurate results than have been obtained here through assumption 
of a hemispherical ground plane. 

Second, since the current distribution measured shows 
deviation from the original sinusoidal form, the measured current 
may be taken as a given current distribution and may be decomposed 
into Fourier Series to add an effect of eachterm. In this way, the 
experimental result may be made to match better against theoretical 
results. 

Third, because the ultimate application of this analysis 
is for a vehicular mounted antenna where a body of a vehicle can be 
taken as an equivalent ground plane of a finite diameter, it is necessary 
to study further to determine how a pa:ticular size of a vehicular body 
and shape and location of the antenna correspond with an antenna ona 
finite disc or semispherical ground plane above an infinite natural 


gound. 


Finally, a detailed study of the measured input impedance 
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given in Chapter 2 and deriving equivaient circuits for each case may 


be useful for design of tuning networks for vehicular mounted antennas. 


APPENDIX A 


RECIPROCITY THEOREM 


Let a source or a distribution of sources maintain a “urrent 


density Jy with associated electromagnetic fields Ey and H, . A 


second source maintains a current density Jo with associated fields 


E, and Hp . Then, according to Lorentz's reciprocal theorem: 


xH,) = J,:E -J,°E 


V- (E, xHy)-V-° 975 1 


2 (A. 1) 


Let's assume that the two sets of fields are defined in free 
Space and are bounded by the same geometrical surfaces. Applying 
volume integration on both sides of Eq. A.1 with dV as volume element, 
the integration performed over a volume V in empty space that is 


bounded by a closed surface S can be written as 
HE Ee ey ay (A. 2) 


where 


SJfG,° B,-5,-B,) av = 0 (A. 3) 


because all volume densities of current vanish in empty space. 


From the divergence theorem 
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ie (E, x H) - (Ex H,) av = g (E, x H, - E, x H,)° dS = 0 
(A. 4) 


Where dS = dS and ft is an external unit normal vector 
of the surface toward volume V. The only resirictions of the two sets 
of fields (E, ; HY) and (E, ; H,) are thit they satisfy Maxwell's equations, 
the volume integration is performed in empty space V and the boundary 
surfaces are in the same geometrical shape for bot sets, but not neces- 


sarily the same physical properties. 


Pe ee ee a Pee ee es 


ak: eae 
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APPENDIX B 


PROOF AND DERIVATION OF EQ. 5. 56 


Equation 20 was written as 


E(x. y, z) = a i, [jwu(n x H) 9 + (nx E)x V'¢+(n- EB) V' olas' 


(5. 50) 


At far-zone, 


r =P (x-x')? + (y-y')? +(z- z')? = IR-R'! (5.29) 


can be simplified as 


r= R- RR = R-R: RF 
Therefore, 
gikR - jkR- Rr 
ane Vaan (5. 55) 
when 
IR| >>IR'l 
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Then 


t 


Nn = i a 
v's a €& v' -jkR -R’ elkR -jKR -R 


: : A 
q R e = -jk R R 


(5. 47) 


’ where the prime notation on operator V indicates that V is operated 


on the primed coordinate system. 


In Eq. 5.50, the radial component of E (x, y, z) canbe 


written as 
jkR ‘KR R' 
4 1 en = = i ae 
; Eg(%) ¥2) = - ge Pe jwu(®x H) - jk(®. nfl : z ds' 
F (B. 1) 


Consider the following integral 


Pay Rr t 
Gg v' (8 ee ). dS = 0 (B.2) 
s 


Where S isa closed surface bounding a volume V. 


Since 


Vx(ga) = Vodxat+gGVxa 


Therefore, Eq. B.2 becomes 


aA = aA 
A eae i. . -1 Py t 
G[S- (x Heri OR +H (Vie ikR- R x H)]as = 0 
Ss 
(B. 4) 
where 
dS = fds 
Knowing 
1A a A = 
vie JKR OR ~jk e JKR 'R R 
and 


Equation B. 4 can be written as 


t 


Pe 
A 'y FW) avd : jKR +R’ 
§ [8 (0 x H eHKR-R (Px H) jie eT IHR RY 


(B. 5) 
Also, from Maxwell's equation 


VixH = -jweg E 


in a source free region. 


Thus, 


Aw =! ~Wr = 
') o E) (~jwe) e JKR “R +R - (Px H) jk @ JiR R] ds' 


fi A =! 
42 G [ieng - x) - i z e KR’ R as = 0 


It is, therefore, true that 


KR: R’ 


GF [iwuy R- hx A) - iK@- By] ec” ds’ = 0 (B.6) 


From Eq. B.1 and B.6. the radial component of E(x, y, z') is proved 
to be zero at far-zone. 
At far-zone, consequently, 


E(x, y,z) = -_ Sf [iwn @x i) 0 J, as’ (5. 56) 


APPENDIX C 


PROOF OF EQ. 5.92 


The purpose of this appendix is to show the steps involved 
in arriving at Eq. 5.92 from Eqs. 5.90 and 5.91, The same steps were 
also taken to get the expression (Eq. 5. 106). 


From Eq. 5.90 and Eq. 5. 91, 


; Egp = U(k? + an ) Bena PT Pa: ha OP nace ose) Pray KR) 
; aR * a ee 
cR? (2n+1)- K+ p._, (ka) (2n+1)- k+ pi 4! (ka) 
: (5, 90) 
1,4 (R) = ae Sin k \d - R) (5.91) 
where 
q d =ath 
Therefore, 
d Ual d n2 
[| Eop1,,(R)dR = ——____*=— ff sink(d- R)(k? +——)p,_ ,(KR)AR 
R=2 (2n+1)-k- Pia (ka) a aR 

Ual d 22 
- — ae ff sink(d- R)(k? + — )p,, , (sR)dR 
(2n+1):k- Prot (ka) a eR? 
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Let 


d 2 
‘ a 
I, = J sin k(d- R)(k? ey Pp, _ ,(R) dR 


(C. 2) 


b 52 
I, = sin k(d- s 
a= J sinka-ri +25 tary an 


aR? 


Changing a variable kR to Z, 


kR = z .’. kdR = dz 
oR? ez? 


The integral I becomes 


kd as 
I, = f{ sink(d- z)(k? +k? £ 
1 

ka 


dz 
az2 ) Py. i (2) eS 


kd 


kd 2 
= k f Sink(d - z) pn. 1 (2) dz+k i sink(d - z) e 
ka ka 


az? 


py. 12) dz (C.3) 


Let 


kd 
Ly = k J sink(d - Zz) p_ 1) dz (C. 4) 
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and 
f sina 2) 2 
I. =k f sink(d- z) p.,(z) dz = (C.5) 
12 ka az? n-1 
Then 
kd 
14 = k f sink(d - z) p,- 2) dz. 
ka 
Let 


1 (2) =u Sink(d- z) dz = dv 


=— 12) = du cosk(d- z) = v 


By part integration 


dp 
I,, = k]cosk(d - z) p 1 (2) - f cosk(d - z) —s— dz 


11 _ 


Integrating the second term of 4 by part again, we get 


kd ap 
af cosk(a - z) —-———- dz 
ka 2 
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ap__ ,(z) a*p___, (ka) 
= | -sink(d - z) ats f sink(d- Zz) —1 4 


z=ka 


The second integral of the above equation is exactly the same 


as Lo- 
Rewriting iy we obtain 
kd 
Li = k eos Zz) eae 1) + sink(d - 2) Pred al - lio 
z=ka 
(C. 6) 


I) = V4 + Lo = [py (kd) - cos k(d - a) Pre 1{ka) - sin k(d - a) Pr-4 a) 


(C.7) 


and 


I, = KK Ps (kd) - cos k(d - a) Pp, (ka) - sin k(d - a) Pas =] (C.8) 


NR ee ae Ee ee ee ee 
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Finally, 
d Ual p (kd) Pp, _ (ka) 
f Eon 1,4 (8) dR = + eae cos k(d - a) dt i ed 
a Py. 1'(ka) Py. 1' (ka) 
P__ 4 (kd) p (ka) 
ge ames k(d- a) B41 (5. 92) 


Pua (Ka Py + 1 (ka) 


APPENDIX D 


NUMERICAL EVALUATION OF BL 


The functions shown in the evaluation of the coefficients B, 
of an infinite series expressing an induced surface current on a hemi- 
spherical ground plane are the weighted spherical hankel functions of 
the second kind andtheir derivatives. 

These functions p(s) can be expressed in the form of a 


series sch as 


A n 
pix) = xh@o) = gre y n+g.k) + 
k=0 (2) x) 


where 


1 . {n + k)! 
n+>>%) = bpa-ks dD) 


T(z+1) = 2! z: integers 
n(x) =f Hay (x) and j = v-1_ 


Also, from the recurrence relations of Hankel functions, the 
following relationship between the spherical Hankel functions and their 


derivatives can be obtained (Ref. 18). 
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fey = a Pe - 222 nen (D. 2) 
and 
dp, (x) 
a ( h xy) = ex) ex ney (D. 3) 
From Eqs. D. 2 and D.3, the following relationships are 
derived 
d p_ (x) 
t pte) = ee = nye) - x jon (D. 4) 
dp_ (x) 
p(x) = Be = xn Pm - nen (D. 5) 
p(x) = xh “)(x) (D. 6) 
n n : 


Equations D.4, D.5, and D. 6 are used with a table for spheri- 


cal Hankel functions to evaluate the coefficients Bo ‘ 


APP. ,NDIX E 


EVALUATION OF R.. in Eq. 5. 168 
1 


Fron the power series representation of cosine integral 
Ci(z), Eq. 5. 165, the terms associated with cosine integral of a small 


argu.uent can be written as 


lim ci(akia + h) . = y+ dn 2k(a+h)+lim fLn(e) (E.1) 
oe) €~O 
lim Ci(2kae) = y+ n(2ka)+lim &nfe) (E. 2) 
€-0 €-0 
lim Cik(2a+h)e = y+ in(k(2a + ns lim £n(e) (E. 3) 
€+0 €+0 
lim Ci(khe) = y+ £n(kh)+lim ff) (E. 4) 
€-0 €-0 


Substituting Eq. E.1 through Eq. E.4 into appropriate places in(1) through 


R= 60 {fn 2- lim ne +y+lim £n(e) + £n(kh) - 7 C08 2k(a + h) £n(2ka) 
r 
1 ame) €-O 
: ; cos 2k(a +h) £n 2k(a + h) + cos 2k(a +h) £n(2a + h)k - Ci(2kh) 


+ = cos 2k(a + h) Ci(4ka) +$ sin 2k(a + h) Si(4ka) - (1 - oe ) sin 2 kh 
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- cos 2k(a + h) cilta + h) - Sin 2k(a + h) Si(2x(20 + m)) 


+4 cos 2k(a + h) Ci (sia + n) + ; sin 2k(a + h) Si (sia + 1) 


(E. 5) 


It should be noticed that the coefficients of y and £n(e) 
associated with the cosine integrals add up to 1. This can be shown 
easily after sO"e algebraic manipulation with trigonometric identities. 


After cancelling out the logarithmic singularities in Eq. E.5 


and collecting terms, the final form of Eq. E.5 can be shown as Eq. 5. 168. 
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11. 
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